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1.0  INTRODUCTION 

1 . 1  General 

The  prediction  and  understanding  of  tide  induced  currents  in  coastal 
waters  and  their  interaction  with  beaches  and  structures  represent  an  intriguing 
and  complex  area  of  research  today.  These  currents  are  especially  significant 
in  the  vicinity  of  tidal  inlets  and  estuary  or  river  mouths.  They  affect 
commerce,  navigation,  small  craft  operations,  water  quality,  bathymetric 
changes  and  shoreline  stability  in  these  densely  populated  coastal  areas. 

From  an  engineering  point  of  view,  the  primary  areas  of  interest  can  be  identi¬ 
fied  as  the  exchange  between  the  interior  and  exterior  waters  connected  by 
tidal  inlets,  the  transport  of  coastal  sediments  and  the  circulation  patterns 
established  in  the  adjacent  waters. 

A  tidal  inlet  is  a  narrow  entrance  region  connecting  an  interior  water 
body  such  as  a  bay  or  lagoon  to  the  exterior  ocean  domain.  Therefore,  inlets 
serve  as  interfaces  between  these  water  bodies  and  contribute  to  their  flushing 
and  exchange  mechanisms,  as  shown  by  Dean  and  Taylor  (1).  They  have  found  that 
the  water  issuing  from  inlets  during  the  ebb  tide  encounters  a  renewal  of  about 
45%  when  it  returns  during  the  following  flood.  This  efficient  mechanism 
of  flushing  is  due  to  the  unique  features  of  the  flow  near  the  narrow  entrance. 
Since  the  actual  hydromechanics  associated  with  entrances  have  not  been  thoroughly 
understood,  prior  modeling  efforts  on  the  circulation  and  mixing  in  bays  have 
often  made  use  of  assumed  boundary  conditions  at  inlets. 

The  mass  transport  near  narrow  entrances  also  has  important  consequences 
with  respect  to  the  transport  of  coastal  sediments.  It  has  long  been  recognized 
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that  tidal  inlets  act  as  sediment  sinks,  trapping  large  amounts  of  littoral 
materials  in  the  form  of  extensive  shoals  encircling  entrance  regions.  There 
is  a  direct  relationship  between  the  deposition/erosion  patterns  near  an  inlet 
and  the  mass  transport  phenomena  associated  with  tidal  currents.  Yet,  the  prior 

|  ] 

design  of  jetties,  navigation  channels,  and  sand  bypass  systems  have  often 
been  based  on  considerations  of  the  poorly  understood  surf  zone  wave  dynamics, 
whereas  the  role  of  tidal  hydromechanics  has  been  largely  ignored.  The  importance 
of  the  latter  effects  has  been  demonstrated  by  Dean  and  Walton  (2).  According 
to  their  study,  the  outer  shoal  volumes  of  major  Florida  inlets  range  between 
106m3  to  2  x  10  0m3,  and  in  some  areas  the  sand  extracted  from  the  adjacent 
beaches  has  caused  40  feet  retrogressions  of  the  shoreline  each  year  for  a 
period  of  50  years. 

The  domain  of  influence  of  the  above  flow  and  mass  transport  phenomena 
extends  from  the  confined  areas  of  bay  or  estuary  waters  to  the  waters  on  the 
semi-infinite  continental  shelf  region.  In  fact,  due  to  the  presence  of  a 
narrow  entrance,  large  scale  circulation  patterns  are  generated  on  both  the 
interior  and  the  exterior  regions  which  are  extremely  important  from  an  oceano¬ 
graphic  point  of  view.  They  affect  the  water  mass  distributions  (Csanady(3)) 
and  mean  circulation  of  adjacent  shallow  seas  (Rouse  and  Coleman  (4)). 

The  preceding  discussion  is  best  understood  when  related  to  tidal 
jets  observed  in  nature.  The  photograph  in  Fig.  1.1(a)  shows  a  starting  tidal 
jet  at  Fort  Pierce  inlet  in  Florida.  The  jet  displays  two  distinct  zones,  a 
frontal  region  characterized  by  vertical  flow  and  a  region  with  primarily 
unidirectional  flow.  The  importance  of  circulation  patterns  induced  by  tidal 
jets  in  adjacent  waters  is  demonstrated  in  Fig.  1.1(b).  This  figure  shows 


n 

Li 


2 


ebbing  flow  at  St.  Mary's  Entrance,  Florida.  The  entrainment  of  effluent 
material  released  far  away  from  the  south  jetty  into  the  jet  is  very  clearly 
seen.  Figure  1 . 2  provides  a  striking  picture  of  two  jets  issuing  from  two  neigh¬ 
boring  inlets  at  Big  Hickory  Pass  and  New  Pass  in  Florida.  The  jets  are  attrac¬ 
ted  and  eventually  attached  to  each  other  due  to  entrainment  velocities.  The 
effects  of  entrainment  into  both  jets  has  caused  recession  in  the  middle  portion 
of  the  beach,  and  hence  offset  the  geometry  of  the  shoreline.  The  smaller  of 
these  two  inlets  has  a  history  of  closing  and  reopening,  and  has  migrated  toward 
the  large  inlet  at  a  rate  of  100  m.  in  6  months. 

The  research  presented  herein  is  concerned  with  the  study  of  circulation 
patterns  induced  on  the  continental  shelf  by  tidal  jets  during  el*b.  Since  they 
occur  in  an  unconfined  environment,  these  flows  are  easier  to  analyze  than  interior 
flows  during  flood  which  occur  in  a  confined  environment.  The  latter  situation 
is  more  complex  and  will  be  treated  in  future  investigations  along  with  the 
mass  transport  phenomena  associated  with  the  flooding  and  ebbing  flows.  Before 
considering  details  of  the  scope  of  present  work,  which  are  presented  in  Sec. 

1.3,  it  is  useful  to  survey  previous  literature  related  to  turbulent  jets  in 
general.  This  is  done  in  the  next  section. 

1.2  Review  of  Pertinent  Literature 

A  bay  or  a  lagoon  region  connected  to  the  ocean  through  a  tidal  inlet  can 
be  Idealized  as  in  Figure  1.3.  The  flow  patterns  associated  with  tidal  currents 
have  rather  unique  features  in  the  vicinity  of  a  narrow  entrance.  As  the  water 
starts  ebbing  into  the  ocean,  an  unsteady  flow  development  is  observed,  as 
schematized  in  Figure  1.3.  Initially,  a  radially  expanding  source  flow  is 
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Definition  Sketch  (a)  Bay-Inlet-Ocean  System  (b)  Surface 
and  Bottom  Variations. 


issued  from  Che  inlet  (Fig.  1.4a).  As  Che  inlet  velocity  increases,  the  flow 
is  separated  from  coastal  boundaries  (Fig.  1.4b)  and  becomes  transformed  from 
a  source  flow  into  a  jet  (Fig.  1.4c).  A  vortex  pair,  shed  off  during  the  im¬ 
pulsive  start  of  the  motion,  moves  offshore,  forming  the  frontal  portion  of 
the  jet.  These  features  of  the  starting  jet  are  shown  photographically  in 
Fig.  1.1a  for  Fort  Pierce  Inlet.  The  structure  of  the  jet  is  determined  by 

the  Reynolds'  number  Re  *  \i\lh  b  /v  where  u  ,  h  ,  b  are  respectively  the 

o»oo  o  o  o 

mean  velocity,  the  depth  and  the  half-width  of  the  inlet,  and  v  the  kinematic 
viscosity.  According  to  Pearce's  (6)  results,  the  flow  is  laminar  for  Re<500, 
becomes  turbulent  in  the  range  Re*1500-2500,  and  remains  fully  turbulent  for 
Re>3000.  Thus,  some  time  after  the  beginning  of  the  ebb  period,  a  fully  developed 
turbulent  jet  is  formed  behind  the  frontal  region. 

The  unsteady  development  of  a  starting  jet  can  approximately  be  represented 
as  a  composite  of  a  steady  turbulent  jet  and  an  unsteady  frontal  region.  This 
assumption  has  been  used  by  Turner  (7),  Tsang  and  Wood  (8),  Tsang  (9),  and 
Middleton  (10)  in  analyzing  buoyant  plumes,  where  plume  solutions  have  been 
coupled  with  thermal  solutions  at  the  frontal  region.  The  thermal  comprises 
a  vortex  ring  in  axisymmetric  starting  plumes,  or  a  vortex  pair  (such  as  shown 
in  Figs.  1.4a  and  1.4b)  in  the  two  dimensional  case.  In  the  case  of  tidal 
jets,  an  order  of  magnitude  analysis  by  Ozsoy  (11)  has  shown  that  unsteady 
effects  are  indeed  negligible  in  the  near  field  of  the  inlet  mouth.  It  has 
further  been  demonstrated  that  the  front  speed  calculations  based  on  steady 
jet  solutions  yield  reasonable  agreement  with  laboratory  observations  of  starting 
jets.  Therefore,  steady  analyses  are  sufficient  to  infer  the  quasi-steady  jet 
properties. 
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A  detailed  review  of  turbulent  jets  can  be  found  in  Albertson  et  al 
(12),  Abramovich  (13),  Schlichting  (14),  Stolzenbach  and  Harleman  (15)  and 
Rajaratnam  (16). 

The  significance  of  jet  flows  at  tidal  inlets  lies  not  only  in  the  dynamics 
and  mass  transport  within  the  jet  itself,  but  also  in  its  effects  on  the  larger 
scale  circulation  patterns.  The  surrounding  fluid  particles  near  the  jet 
boundaries  are  entrained  and  carried  away  by  the  jet,  owing  to  the  friction 
between  the  two  zones  (13,  14).  The  resulting  velocity  of  entrainment  causes 
induced  circulation  in  the  outer  part  of  the  fluid  away  from  the  jet  region  (17). 
Such  entrainment-induced  flows  are  possible  causes  for  gyres  generated  in  the 
vicinity  of  river  mouths  as  reported  by  Wiseman  et  al  (5)  and  circulations 
observed  by  Olsen  (18)  near  tidal  inlets. 

A  first  attempt  has  been  made  by  French  (19)  to  study  the  ocean  flow 
patterns  during  both  portions  of  the  tidal  cycle.  The  flood  flow  is  analyzed 
by  making  use  of  the  free  streamline  theory,  including  the  effects  of  long 
jetties.  The  ebb  flow  is  represented  by  Tollmien's  solution  for  two-dimensional 
turbulent  jets,  neglecting  the  effects  of  bottom  friction  and  variable  bathy¬ 
metry.  In  reality,  both  are  significant  factors  in  the  analysis  of  jets. 

The  effects  of  bottom  friction  and  depth  variations  on  a  jet  have  been 
first  taken  into  account  by  Borichansky  and  Mikhailov  (20),  however,  their 
results  are  erroneous,  since  the  jet  momentum  has  been  assumed  to  be  conserved, 
even  in  the  presence  of  bottom  friction.  In  spite  of  these  and  other  deficiencies 
in  the  analysis,  these  writers  have  shown  the  exponential  dependence  of  jet 
expansion  and  velocity  decay  on  offshore  distance.  Furthermore,  they  have  demon¬ 
strated  the  opposing  effects  of  a  seaward  slope  and  the  bottom  friction. 
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A  similar  analysis,  including  only  the  bottom  frictional  effects,  has 
been  performed  by  Taylor  and  Dean  (21).  Upon  neglecting  lateral  mixing  and 
entrainment,  they  have  obtained  solutions  for  the  jet  half-width  b  and  the  jet 
velocity  u  (assumed  uniform  across  the  jet)  as 


iL.  ,  .«/ 2 

b 


(1.1a) 


u 

u 


-Ui/2 


(1.1b) 


where  p  =  fb  /8h  is  the  parameter  specifying  bottom  friction  effects  (f  being 
o  o 

a  friction  coefficient),  and  £  *  x/b  the  normalized  offshore  distance  from 

o 

the  inlet.  Note  that  the  aspect  ratio  bQ/hQ  is  also  incorporated  into  the 
parameter  u,  such  that  increasing  the  ratio  of  inlet  width  to  depth,  in  effect, 
increases  the  frictional  resistance. 

A  more  detailed  analysis  of  the  jet  has  been  performed  by  Unluata  and 
Ozsoy  (22)  and  Ozsoy  (11),  in  which  general  solutions  have  been  obtained  in¬ 
cluding  the  effects  of  bottom  friction,  bathymetric  changes,  turbulent  mixing 
and  lateral  entrainment.  The  velocity  distributions  are  assumed  to  be  self¬ 
similar  and  the  boundary  layer  approximations  are  utilized  in  obtaining  analytical 
solutions  to  the  governing  Navier-Stok.es  equations.  In  the  near  field  of  the 
inlet.  It  is  shown  that  the  unsteady  effects,  the  Coriolis  forces,  and  the  surface 
gradients  are  negligible.  The  velocity  vg  of  entrainment  of  fluid  particles 

into  the  jet  is  assumed  to  be  v  «  ciu  after  Morton  et  al  (23),  where  a  is 

c  c 

the  entrainment  coefficient  (o  »  0.05),  and  u  the  jet  centerline  velocity. 

c 

The  analysis  is,  in  general,  similar  to  the  one  used  by  Stolzenbach  and  Harleman 
(15).  In  the  absence  of  depth  variations,  and  for  large  distance  (£-*«>),  these 


solutions  are  simplified  to  yield 
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where  k^,  are  constants.  The  general  features  of  these  solutions  are 
shown  in  Fig.  3.2  of  Sec.  3.2  which  contains  details  of  Ozsoy's  (11)  work. 

The  exponential  dependence  of  jet  width  and  velocity  on  distance  here  is 
considerably  different  from  that  obtained  by  Taylor  and  Dean  (21)  in  Eqs.  (1.1), 
since  the  lateral  mixing  and  entrainment  have  been  neglected  by  the  latter 
authors.  Due  to  the  additional  effects  of  lateral  entrainment,  the  expansion  and 
velocity  decay  rates  are  doubled  in  Eqs.  (1.2)  as  compared  to  the  bottom 
frictional  jet  solution  of  Ref.  (21).  It  is  also  very  important  that  the  exponen¬ 
tial  behaviour  of  the  bottom  frictional  jet  differs  significantly  from  the 
linear  expansion  (b~C)  of  a  classical  two  dimensional  jet  (12,  13,  14),  in  which 
the  centerline  velocity  decays  as  uc~5  •  In  the  absence  of  bottom  friction 
(p»0),  the  general  solutions  of  Unluata  and  Ozsoy  (22)  reduce  identically  to 
the  classical  jet  solutions.  However,  even  the  linear  jet  growth  in  the  classical 
theory  has  been  questioned  by  Kotsovinos  (24)  who  showed  that  in  most  experiments 
a  nonlinear  growth  is  observed.  Kotsovinos  (24)  attributed  this  nonlinear  growth 
to  ambient  turbulence  levels,  but  the  frictional  effects  of  side  plates  used 
in  the  experiments  seem  also  plausible,  based  on  the  above  arguments.  In  his  further 
research,  Kotsovinos  (25)  has  demonstrated  that  the  constancy  of  momentum  flux  in 
the  axial  direction  of  the  jet  is  an  approximation  (even  in  the  absence  of  bottom 
friction)  because  the  jet-induced  flow  has  a  component  opposite  to  the  main 
jet  flow  and  because  of  the  pressure  field  generated  in  the  ambient  fluid. 


Unluata  and  Ozsoy  (22)  and  Ozsoy  (11)  have  also  shown  the  counteracting 
effects  of  a  seaward  slope  and  the  bottom  friction  In  determining  the  jet 
behaviour.  In  particular,  when  these  opposing  effects  are  in  balance,  a  linear 
growth  similar  to  the  classical  jet  is  recovered  due  to  entrainment  alone. 

For  bottom  slopes  in  excess  of  this  balanced  value,  the  growth  becomes  logarith¬ 
mic  representing  a  lower  expansion  rate.  The  effects  of  bottom  slope  on  the  jet 

development  is  shown  in  Fig.  1.5,  where  u  •  mb  /h  is  a  slope  parameter,  m 

o  o 

being  the  bottom  slope.  In  the  case  of  arbitrary  depth  variations:  h*h(x),  the 
jet  is  found  to  go  through  a  series  of  expansions  and  contractions,  depending 
on  the  local  variations  in  the  bathymetry.  One  such  example  is  Jupiter 
Inlet,  Florida  shown  in  Fig.  1.6.  Analytical  solutions  obtained  for  this  case 
are  compared  with  observations  in  Fig.  1.7.  The  jet  expands  over  the  shoals 
and  contracts  afterwards  due  to  depth  increases. 

The  effects  of  cross-currents  in  the  ambient  waters  have  also  been  included 
by  Ozsoy  (11).  The  centerline  loci  and  the  half-width  for  these  deflected  jets 
are  shown  in  Fig.  1.8,  where  U.  *  u  /u  ,  u  being  the  cross-current  velocity. 

A  <1  O  cl 

Examples  of  deflected  jets  are  provided  by  South  Lake  Worth  and  Baker's 
Haulover  Inlets  in  Florida. 

The  offshore  extent  of  tidal  jets  is  often  in  the  order  of  L  -  10  km. 

u 

A  simple  calculation  of  the  Rossby  number  yields  R  =  *  0(1),  indicating  the 

o  pL 

potential  relevance  of  rotational  effects.  The  effects  of  earth's  rotation  on 
horizontal  jets  have  been  investigated  by  Gadgil  (26)  and  Savage  and  Sobey  (27). 
In  the  case  of  shallow  water  jets  (jets  covering  the  full  depth),  analytical  and 
experimental  studies  by  Savage  and  Sobey  (26)  indicate  that  the  jet  path  remains 
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Figure  1.7  Prototype  Calculations  for  the  Jet  at  Jupiter  Inlet,  Florida 

(a)  Depth  Variations  (b)  Jet  Half-Width  (c)  Centerline  Velocity 
(Courtesy  Ref.  (11)) 
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straight  in  spite  of  the  Coriolis  forces.  This  appears  to  be  a  confirmation  of 
the  Taylor-Proudman  theorem  (Greenspan  (28))  in  the  theory  of  rotating  fluids, 
since  the  vertical  velocity  distributions  are  not  affected  by  rotation.  However, 
the  rotational  effects  are  felt  through  the  development  of  Ekman  boundary  layers, 
which  effectively  increases  bottom  friction.  These  frictional  losses  induced 
by  rotation  have  been  analyzed  by  Gadgil  (26),  for  laminar,  quasi-geostrophic 
jets.  Her  results  indicate  that  there  exists  a  finite  penetration  distance  in 
which  all  jet  momentum  is  lost.  As  opposed  to  shallow  water  (or  two  dimen¬ 
sional)  horizontal  jets,  jets  in  deep  water  have  completely  different  characteris¬ 
tics  (27),  since  they  are  not  able  to  support  rotational  pressure  gradients 
and  the  jet  path  curves  in  a  spiral  shape.  Such  deep  water  jets  can  be  formed 
when  the  orifice  dimensions  are  much  smaller  than  the  depth  or,  in  the  case  of 
tidal  effluents,  when  the  flow  separates  from  the  bottom  due  to  either  buoyancy 
or  a  steep  slope.  The  combined  effects  of  unsteady  growth  and  rotation  may 
also  deflect  the  jet  during  its  initial  stages  as  shown  by  Savage  and  Sobey  (27) . 
In  the  presence  of  depth  variations,  rotational  effects  may  also  become  quite 
significant  due  to  the  conservation  of  potential  vorticity  as  shown  in  (28). 

Due  to  fresh  water  discharges  into  bay  or  lagoon  waters  adjacent  to 
tidal  inlets,  tidal  effluents  can  often  be  buoyant.  In  fact,  such  jet  flows 
in  the  presence  of  stratification  have  been  studied  in  the  field  and  documented 
by  Wright  et  al  (29)  and  Wright  and  Sonu  (30)  at  tidal  inlets  and  by  Wright 
and  Coleman  (31,  32),  and  Garvine  (33)  at  river  mouths.  Due  to  buoyancy,  the 
plumes  separate  from  the  bottom,  and  therefore  bottom  friction  and  bathymetric 
variations  lose  their  importance  at  some  distance  offshore  from  the  inlet. 

However,  in  the  near  field  of  the  mouth,  bottom  topography  and  friction  play 


significant  roles,  as  shown  in  (32,  30,  34).  In  the  presence  of  density  con¬ 
trasts,  the  most  important  mechanism  in  the  far  field  is  probably  the  buoyant 
spreading  of  the  plume  at  the  expense  of  a  decrease  in  the  effluent  depth. 

By  neglecting  the  entrainment  and  mixing  processes,  Bondar  (35)  has  derived 
expressions  for  the  plume  half-width  b  and  interfacial  depth  h: 

=  (1+*C)2/3  (1.3a) 

U 

O 

~r  -  (1+<K)'2/3  (l-3b) 

o 

where  $  (2-y)/u  is  a  measure  of  the  ratio  of  the  buoyant  to  inertial 

o  o 

forces  (proportional  to  the  densimetric  Froude  number),  and  y  =  1  -  Pj/Ps> 
a  parameter  comparing  the  fresh  and  salt  water  densities  and  ps,  respectively. 

The  plume  velocity  u  is  constant  in  this  model  (u  =  u^)  since  the  turbulent  mixing 
with  the  ambient  waters  is  neglected.  Observations  by  Wright  and  Coleman  (32) 
at  the  mouth  of  South  Pass,  Mississippi  River  Delta  support  the  predictions 
based  on  Eqs.  (1.3).  Calculations  by  Wright  and  Coleman  (31)  inc’uding  vertical 
entrainment  and  mixing  at  the  interface  also  indicate  an  approximate  agreement 
with  observations.  Takano  (35)  has  neglected  the  inertial  effects,  but  has 
included  buoyancy  and  shear  forces  as  well  as  entrainment  in  his  analysis,  and 
has  obtained 
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where  h  is  the  interfacial  depth  at  the  plume  centerline.  A  complete  analysis 
including  all  the  near  field  processes  such  as  mixing,  entrainment  and  buoyant 
spreading  has  been  performed  by  Stolzenbach  and  Harleman  (15).  Comparison  of 
these  results  with  laboratory  experiments  show  that  the  model  is  sufficiently 
accurate  in  predicting  both  near  and  far  field  properties  of  buoyant  surface 
plumes  in  the  absence  of  bottom  interaction. 

The  effects  of  cross  currents  on  buoyant  plumes  have  been  analyzed  by 
Stolzenbach  and  Harleman  (15).  The  effects  of  earth's  rotation  have  been  taken 
into  account  by  Takano  (35).  Numerical  modeling  by  Waldrop  and  Farmer  (36) 
has  revealed  the  formation  of  vertical  secondary  circulations  due  to  buoyant 
spreading.  Unsteady  tidal  effects  on  river  plumes  have  been  discussed  by 
Wright  and  Coleman  (32)  and  studied  by  Yoshida  and  Kashiwamura  (37). 

The  interaction  between  incident  waves  and  river  plumes  has  been  studied  by 
Wright  (34).  Due  to  wave  transformation  and  breaking,  a  rapid  deceleration, 
vertical  mixing  and  abrupt  widening  of  the  plume  is  observed,  and  a  correlation 
has  been  obtained  between  the  significant  wave  height  variations  and  the  river 
stage.  However,  the  latter  argument  is  weakened  due  to  the  presence  of  river 
flood  conditions  at  the  time  of  the  experiments. 

The  unsteady  development  of  a  buoyant  plume  issuing  from  a  tidal  inlet  has 
been  described  by  Wright  and  Sonu  (30),  and  Sonu  and  Wright  (34).  The  latter 
study  has  also  revealed  a  cross-current  bypassing  the  plume  under  it. 

The  outer  circulation  generated  by  a  tidal  jet  or  a  river  plume  has  received 
little  attention  to  date,  even  though  it  is  one  of  the  basic  characteristics 
of  jets  or  plumes.  The  streamlines  for  the  entrainment  induced  outer  flow  can 
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be  analyzed  by  using  potential  flow  theory  in  the  domain  external  to  the  jet 
and  by  replacing  the  Jet  with  a  line  sink  distribution;  however,  such  analyses 
have  been  largely  ignored  in  the  previous  literature  on  tidal  or  river  plumes. 

The  above  method  has  been  first  introduced  by  Taylor  (38)  (see  Turner  (17))  to 
analyze  the  outer  flow  due  to  a  buoyant  plume.  It  is  also  conceivable  in  the 
case  of  unsteady  tidal  jets  shown  in  Figures  1.2a  and  1 . 2b  that  the  propagation 
of  the  jet  across  the  continental  shelf  region  sets  up  large  scale  unsteady 
circulation  patterns.  Similar  outer  circulation  patterns  due  to  a  plane 
starting  plume  have  been  studied  by  Tsang  and  Wood  (8),  by  replacing  the  main 
part  of  the  plume  by  a  line  sink  and  superposing  a  point  source  and  a  vortex 
doublet  at  the  front.  The  general  flow  patterns  are  similar  to  the  case  of  a 
tidal  jet,  but  no  literature  exists  either  to  document  or  study  these  effects. 

1 . 3  Scope  of  Present  Investigation 

The  objective  of  this  Investigation  was  to  studv  the  large  scale  circulations 
on  the  continental  shelf  region  induced  by  tidal  Jets. 

The  results  of  research  presented  herein  are  concerned  mainlv  with  steady 
jets  issuing  into  a  shallow,  vertically  well-mixed  bodv  of  water.  The 
Jets  are  assumed  to  be  isothermal  and  salinity  gradients  are  assumed  to  be 
absent  so  that  buoyancy  effects  are  excluded.  In  this  initial  phase  of 
9tudv,  the  considerations  of  unsteadiness,  the  Earth's  rotation,  wave 
climate,  and  a  confined  environment  are  ignored;  however,  the  important 
effects  of  bottom  friction  in  the  jet  have  been  Included.  Ozsoy's  (11)  analysis 
has  shown  that  in  the  near  field  of  the  Jet  Coriolis  forces  are  negligible,  and 
relative  to  the  time  scale  of  interest  (30-60  mins,  typically),  the  unsteady 
effects  are  unimportant.  Furthermore,  factors  such  as  buoyancy,  wave  processes. 
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and  pressure  gradients  (due  to  a  jet  issuing  into  the  confined  environment  of  a 
bay  during  a  flood  tide)  constitute  additional  complexities  and  are  not  within 
the  scope  of  the  work  presented  here.  Finally,  during  major  portions  of  the  study, 
the  jet  has  been  assumed  to  discharge  over  constant  bathymetry,  as  this  results 
in  significant  analytical  simplifications  and  also  because  a  meaningful  first 
insight  into  the  solution  behaviour  is  obtained.  This  assumption  was  relaxed 
subsequently  and  some  preliminary  elementary  solutions,  which  are  building  blocks 
for  determining  Jet-induced  circulation,  have  been  derived. 

The  mathematical  formulation  of  circulation  generated  by  a  jet  issuing 
into  a  shallow  ebbing  ocean  is  described  in  Sec.  2.0.  A  general  partial 
differential  equation  governing  steady  and  irrotational  motion  of  an  inviscid 
fluid  over  arbitrary  bottom  topography  has  been  derived  first.  When  specialized 
to  the  case  of  constant  bottom,  this  equation  reduces  to  Laplace's  equation  for 
potential  flow.  Since  Laplace's  equation  possesses  well-known  solutions,  such 
as  sink  flows,  the  jet  which  entrains  the  surrounding  fluid  was  replaced  by  an 
equivalent  sink  distribution  alongs  its  axis.  This  approach  is  due  to  Taylor 
(11).  The  details  of  this  procedure  are  contained  in  Sec.  3.0  and  the  flow 
field  induced  by  a  simple  jet  (perpendicular  to  the  coastline)  has  also  been 
solved  in  that  section.  These  analyses  can  be  readily  extended  to  treat 
other  flow  situations  related  to  tidal  inlets  which  occur  in  nature.  For 
example,  a  tidal  jet  discharging  through  a  jetty,  two  adjacent  jets  as  in  Fig. 

1.2,  and  a  jet  which  is  not  perpendicular  to  the  coastline.  Sec.  4.0  describes  the 
solutions  for  these  cases.  Numerical  results  for  circulation  induced  by  jets 
are  given  in  Sec.  5.0  in  terms  of  streamline  patterns  and  variations  of  the 
alongshore  current.  Formulation  of  the  variable  bathymetry  problem  is  presented 
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in  Sec.  6.0.  An  elementary  singularity  solution  equivalent  to  the  sink  for 
Laplace's  equation,  has  been  obtained  for  the  special  case  of  an  exponentially 
sloping  bottom.  Application  of  this  solution  to  the  determination  of  jet-induced 
flow  has  been  outlined.  Sec.  6.0  also  contains  a  preliminary  formulation  of  the 
unsteady  circulation  problem  which  takes  into  account  the  effects  of  Earth's 
rotation,  and  bathymetry  changes.  Finally,  in  Sec.  7.0,  conclusions  derived  from 
the  present  research  effort  have  been  described  and  some  suggestions  for  future 


studies  in  this  area  have  been  advanced. 


2.0  MATHEMATICAL  FORMULATION  OF  JET- INDUCED  CIRCULATION 


2. I  Depth-  and  Time- averaged  Governing  Equations 


The  fluid  motion  in  the  bay/ocean  system  shown  in  Figure  1.3a  is  described 
in  the  most  general  sense  by  the  well-known  unsteady  incompressible  Navier-Stokes 
equations  in  three  dimensions.  The  waters  in  the  vicinity  of  tidal  inlets  are 
generally  shallow  and  the  time  scale  of  unsteadiness  equals  the  period  (T)  of 
tidal  variations.  Therefore,  it  is  possible  to  simplify  the  Navier-Stokes  equations 
by  (i)  averaging  over  the  depth  (z  direction)  and  (ii)  averaging  over  a  time 
period  much  longer  compared  to  the  time  scale  of  turbulent  fluctuations  but  much 
shorter  compared  to  the  period  of  tidal  variations.  The  details  of  these  operations 
are  given  in  Ref.  (11)  and  the  resulting  two  dimensional  (x  and  y)  turbulent  motion 
is  governed  by  the  following  equations. 
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where  n  is  the  free  surface  elevation  and  H  *  n  +  h  is  the  total  depth  (Figure 
1.3b).  The  depth-averaged  horizontal  velocity  components  are  denoted  by  u  and  v  and 
the  magnitude  of  velocity  is  given  by  q  *  (u2  +  v2)^.  In  Eqs.  (2.2)  and  (2.3),8,  p, 
f ,  and  g  represent  the  Coriolis  parameter,  density,  bottom  friction  factor,  and 
gravitational  acceleration  respectively.  Fxx»  Fxv  "  ^vX*  anc*  Fyv  are  comPonents 
of  the  depth-integrated  stress  tensor  (including  Reynolds  stresses). 
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Eqs.  (2.1),  (2.2),  and  (2.3)  present  a  formidable  mathematical  problem  which 
cannot  be  solved  analytically  to  obtain  a  general  solution  in  closed  form.  These 
equations  are  simultaneous  partial  differential  equations  in  3  unknowns  u,  v,  and  n, 
and  must  be  solved  numerically  on  a  digital  computer  using  finite  difference 
or  finite  element  techniques.  Such  an  effort,  however,  constitutes  a  major  under¬ 
taking  and  is  not  justified  for  the  fundamental  research  investigation  reported 
herein.  Furthermore,  one  can  easily  get  involved  in  complications  of  numerics 
and  lose  sight  of  the  physics,  thereby  making  it  difficult  to  ascertain  the  correct¬ 
ness  of  the  solutions  obtained.  Lastly,  and  most  importantly,  the  costs  associated 
with  a  complete  numerical  approach  are  very  high.  The  foregoing  considerations 
led  to  the  treatment  of  the  fluid  motion  as  a  steady  potential  flow  and  is  described 
in  the  next  section. 

2. 2  Shelf  Circulation  As  A  Steady  Potential  Motion 

Before  addressing  the  complex  problem  of  shelf  circulation  motions  which 
include  the  effects  of  unsteadiness,  bottom  friction,  and  Earth's  rotation,  it 
is  useful  to  introduce  the  assumption  of  steady  potential  flow.  In  other  words, 
the  motion  induced  by  a  tidal  jet  in  the  surrounding  fluid  is  assumed  to  be  steady, 
inviscid  and  irrotational,  and  in  addition,  the  effects  of  Earth's  rotation  are 
ignored.  It  may  be  pointed  out  that  the  latter  simplification  rules  out,  a 
priori,  the  generation  of  vorticity  (even  in  an  inviscid  fluid)  as  a  result  of 
bathymetric  variations.  The  steady  potential  motion  offers  the  following  advantages 

o  A  simplified  mathematical  formulation  which  preserves  the  essential 
physics  of  the  problem 
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o 


The  steadiness  of  motion  makes  It  possible  to  use  the  existing  theory 
of  a  turbulent  jet  with  bottom  friction  and  bottom  topography  (11). 

Thus,  although  bottom  friction  is  excluded  in  the  analysis  of  the  outer 
flow,  it  is  included  in  the  inner  jet  flow. 

o  Steady  potential  flow  is  governed  by  a  single  linear  partial  differential 
equation  which  is  amenable  to  analytical  treatment.  This  equation 
contains  partial  derivatives  with  variable  coefficients  for  arbitrary 
bottom  topography  and  reduces  to  Laplace's  equation  for  constant 
bathymetry.  A  large  body  of  knowledge  is  then  readily  available  through 
the  well-known  solutions  of  Laplace's  equation. 

o  Solution  of  the  relatively  simple  potential  flow  problem  yields  a  valuable 
insight  which  will  be  useful  in  future  studies  of  more  complex  problems* 
taking  into  account  the  effects  of  unsteadiness,  Earth's  rotation,  and 
bottom  friction. 

It  should  be  brought  to  the  reader's  attention  that  the  term  "steadiness"  in 
the  above  discussion  implies  time-invarient  motion  over  a  period  which  is  large 
compared  to  the  turbulence  time  scale  but  small  compared  to  the  tidal  period 
(e.g.  30  minutes,  typically).  The  derivation  of  equations  governing  steady  poten¬ 
tial  flow  with  variable  bathymetry  is  presented  below. 

Under  steady  flow  conditions  3q/3t  vanishes  in  Eq.  (2.1)  and  it  reduces  to 

KHa)  |i«v) 

3x  3y 
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Assuming  further  that  the  displacement  n  of  the  free  surface  Is  small  compared 


to  the  Stillwater  depth  h  (Figure  1.3)  at  every  point  In  the  flow  field,  l.e. , 
n  <<  h,  Eq.  (2.4)  becomes, 

MhU>^hV)  ■«  <2.5) 


where  the  function  h  (x,y)  represents  bottom  variations.  Eq.  (2.5)  now  contains 
two  unknowns,  u  and  v,  and  it  is  independent  of  the  free  surface  variations  n- 
Therefore,  it  can  be  effectively  decoupled  from  Eqs.  (2.2)  and  (2.3)  if  an 
additional  relation  between  u  and  v  is  available.  Another  equation  in  these 
unknowns  is  provided  by  the  condition  of  irrotationality,  under  the  assumption 
of  potential  motion,  and  it  reads 


3u  3v 
3y  3x 


0 


(2.6) 


for  two  dimensional  flow.  Eq.  (2.6)  implies  the  existence  of  a  velocity  potential 
function  $  (x,y)  such  that 


u 


34> 

3x 


(2.7a) 


v  ■ 


M. 

3y 


(2.7b) 


Substituting  Eq.  (2.7)  into  Eq.  (2.6),  and  simplifying. 


h  h 

+  \T  *x  +  lT  *y  "  0  (2‘8) 

This  equation  governs  the  steady  potential  motion  of  the  continental  shelf  waters. 
It  is  a  linear  partial  differential  equation  of  second  order  with  variable  coef¬ 
ficients  h  /h  and  tWh.  Eq.  (2.8)  brings  out  clearly  the  Laplacean  operator  on 
the  left  hand  side;  and  together  with  the  two  additional  terms  it  represents  a 
generalized  elliptic  equation  in  the  plane. 
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2.3  Steady  Potential  Motion  With  Constant  Bathymetry 


Eq.  (2.8)  suggests  a  very  attractive  simplification,  that  of  constant  or 

flat  bottom,  for  which  h  *  h  30.  The  result  is  Laplace's  equation 

x  y 

72<fr  *  0  (2.9) 

For  the  present  investigation,  solutions  of  Eq.  (2.9)  which  have  been  studied  in 
detail,  must  be  coupled  to  the  characteristics  of  the  tidal  jet.  The  assumption 
of  constant  bathymetry  allows  the  jet  characteristics  (such  as  centerline  velocity 
and  width)  to  be  expressed  in  closed  form.  Clearly,  important  mathematical  simpli¬ 
fications  result  for  the  case  of  constant  bathymetry.  From  the  physical  view¬ 
point,  a  number  of  tidal  inlets  display  very  gradually  varying  bottom  topography 
and  may  be  treated  under  the  above  formulation  as  a  first  approximation. 

The  present  work  is  mainly  concerned  with,  but  not  limited  to,  the  study  of 
continental  shelf  circulation  for  constant  bathymetry.  A  preliminary  formulation 
including  the  effects  of  variable  topography,  unsteadiness,  and  Earth's  rotation 
is  given  in  Section  6.0. 

2.4  Stream  Function  Formulation  Of  Steady  Potential  Motion 

The  problem  formulation  presented  in  the  previous  section  was  in  terms  of  a 
velocity  potential  <p.  An  equivalent  formulation  in  terms  of  a  stream  function  i p 
is  also  possible.  It  may  be  recalled  from  elementary  fluid  mechanics  that  the 
curves  of  constant  $  are  orthogonal  to  curves  of  constant  \ h  everywhere  in  the  flow 
field.  The  basic  kinematical  equations  of  steady  and  irrotational  motion  are,  of 
course,  the  same  as  before,  viz.  Eqs.  (2.5)  and  (2.6).  Instead  of  defining  a 
velocity  potential  through  Eqs.  (2.7),  a  'modified'  stream  function  is  defined 
as  follows. 
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1  3U< 
U  "  "  h  3y 


(2.10a) 


v  -  iii 

h  3x 


(2.10b) 


Notice  that  these  equations  satisfy  the  continuity  Eq.  (2.5)  identically.  Sub¬ 
stituting  the  definitions  in  Eqs.  (2. 10a)  and  (2. 10b)  into  Eq.  (2.6),  and  simplifying, 
the  governing  equation  for  a  steady  potential  motion  is 


h  h 

V2ip - £  ill - ^  1 1>  *0 

h  x  h  y 


(2.11) 


The  similarity  between  Eqs.  (2.11)  and  (2.8)  is  readily  observed  and  the  comments 
following  the  latter  equation  are  equally  applicable. 

For  the  special  case  of  constant  bathymetry,  Eq.  (2.11)  reduces  to  a  Laplace's 
equation  in  iji,  or. 


h  -  0 


(2.12) 


Also,  in  this  case,  h  drops  out  of  the  differentiations  in  Eq.  (2.5),  and  the 
modified  stream  functions  in  Eqs.  (2.10a)  and  (2. 10b)  can  be  replaced  by  the 
standard  definitions. 


u--|t 

3y 


(2.13a) 


(2.13b) 


There  are  now  two  alternate  formulations,  one  in  terms  of  a  velocity  potential 
given  by  Eqs.  (2.9),  (2.7a),  and  (2.7b),  and  the  other  in  terms  of  a  stream  function, 
represented  by  Eqs.  (2.12),  (2.13a),  and  (2.13b).  The  selection  of  one  approach 
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over  the  other  is  a  matter  of  convenience  of  mathematical  analysis  and  inter¬ 
pretation  of  results.  In  the  present  work,  it  was  found  useful  to  employ  the 
stream  function  formulation  to  study  the  constant  bathymetry  problem.  On  the  other 
hand,  a  velocity  potential  formulation  appears  to  be  more  meaningful  and  illumi¬ 
nating  for  the  variable  bathymetry  case  (Section  6.0).  The  constant  bottom 
topography  solutions  are  discussed  first  in  Sec.  3.0.  The  coupling  between 
solutions  of  Eq.  (2.13)  and  the  jet  characteristics  is  also  discussed  in  that 


3.0  STEADY  POTENTIAL  FLOW  WITH  CONSTANT  BATHYMETRY  -  THE  SIMPLE  JET 


3. 1  Representation  of  a  Tidal  Jet  By  a  Sink  Distribution 

It  has  been  shown  in  the  previous  section  that  the  steady  potential  flow 
induced  by  a  tidal  jet  is  governed  by  Laplace's  equation  (2.12)  when  the  bottom 
topography  is  assumed  to  be  constant.  From  potential  theory,  it  is  well-known 
that  Laplace's  equation  can  be  solved  in  terms  of  elementary  solutions  which 
correspond  to  various  physical  flow  situations  (39).  In  the  present  investiga¬ 
tion,  the  jet  emerges  into  a  body  of  water  which  is  governed  by  potential  motion, 
and  therefore,  only  the  jet  need  be  represented  in  terms  of  elementary  solutions. 
Thus,  the  treatment  of  jet-induced  circulation  reduces  to  relating  the  jet  charac¬ 
teristics  to  the  elementary  solutions.  This  idea  is  central  to  the  work  described 
in  this  report  and  was  originally  introduced  by  Taylor  (38)  during  his  study  of 
flow  fields  induced  by  steady  buoyant  plumes.  The  plume  was  replaced  by  a  dis¬ 
tribution  of  sink  singularities  which  constitute  elementary  solutions  of  Eq. 

(2.12)  and  the  external  motion  was  assumed  to  be  potential.  Following  this  appoach, 
the  tidal  jet  was  represented  by  a  distribution  of  sinks  along  its  axis  and  as 
outlined  before,  the  surrounding  frictionless  fluid  motion  over  a  constant  bottom 
was  treated  as  potential.  Thus  the  entire  flow  field  satisfies  Laplaces  equation, 
except  for  the  jet  axis  itself. 

The  representation  of  the  jet  by  a  sink  distribution  also  makes  physical 
sense.  It  is  an  observed  fact  that  a  tidal  jet  draws  in  fluid  adjacent  to 
it  (Fig.  1.1b),  and  thus,  it  effectively  acts  much  like  a  distribution  of  sinks. 

Of  course,  it  is  necessary  to  establish  an  appropriate  correspondence  between  the 
parameters  that  characterize  a  jet  and  those  which  characterize  a  sink  distribution. 
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This  discussion  is  deferred  to  Sec.  3.3.  Before  deriving  the  mathematical  ex¬ 
pressions  for  the  jet-induced  motion,  it  is  important  to  briefly  review  the  jet 
solution.  The  following  section  is  devoted  to  that  purpose. 

3.2  Review  Of  Jet  Solution  For  Constant  Bathymetry 

Ozsoy  (11)  has  obtained  an  analytical  solution  for  a  two  dimensional 
turbulent  jet  with  bottom  friction  and  arbitrary  bathymetry.  Figure  3.1  shows 
the  jet  configuration  considered  in  Ref.  (11).  The  general  expressions  for  the  jet 
characteristics  of  interest;  namely,  core  half-width  and  jet  half-width  in  the  zone 
of  flow  establishment  (ZOFE),  and  centerline  velocity  and  jet  width  in  the  zone  of 
established  flow  (ZOEF).  Whereas  the  solutions  for  arbitrary  bathymetry  must 
be  obtained  numerically,  closed  form  solutions  can  be  derived  for  the  soecial 
case  of  constant  depth.  The  solutions  are  given  below. 

ZOFE,  C  <  5  : 

s 

I,  e'MC  -  I-  (1  +  a?) 

R(5>  -  — — r — r -  <3.0 

-  L2 


B(C) 


<1-I2)  (1+aC)  ~  (1-Ij)  e 


(3.2) 


(3.3) 
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Uo  mwr 


uc(x)=  u(x,o) 


ZOFE 


ZOEF 


Figure  3.1  Definition  Sketch  For  Shallow  Water  Jet  Flows 
(Courtesy  Ref.  (11)) 


where  £,  R,  B,  and  U  are  offshore  distance,  core  half-width,  jet  half-width,  and 
centerline  velocity,  respectively.  These  quantities  are  dimensionless  and  are 
defined  as  follows. 


£  =  x/b 

o 

R(0  -  f 

o 

BU)  -f 

o 

HO 

o 

Furthermore,  u  is  a  modified  friction  parameter  defined  by 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


(3.9) 


with  f  *  Darcy-Weisbach  friction  coefficient.  In  Eqs.  (3.1)  -  (3.4),  1^  -  0.450 
and  I,  -  0.316.  The  core  length  £  Is  given  by  the  condition  R(£)  *  0  in  Eq.  (3.1). 

Eqs.  (3.1)  -  (3.4)  reveal  that,  in  the  presence  of  friction,  the  jet  width  grows 

and  the  centerline  velocity  decays  exponentially  at  sufficiently  large  distances 

away  from  the  inlet  as  opposed  to  the  linear  growth  of  the  classical  two-dimensional 

_u 

turbulent  jet  in  which  the  centerline  velocity  decays  as  £  .  As  with  the 

classical  two  dimensional  jet,  the  linear  behaviour  of  the  jet  width  occurs 
only  near  the  inlet  and  in  the  core  region.  Figure  3.2  illustrates  the  variations 
of  Jet  half-width  and  centerline  velocity  against  distance  from  the  inlet  with  p 
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as  a  parameter.  As  seen  from  Eq.  (3.9),  p  combines  the  relative  effects  of  bottom 

friction  and  inlet  aspect  ratio  b  /h  .  When  friction  is  large  or  the  inlet  width 

o  o 

is  large  compared  to  depth,  the  jet  merely  explodes  as  it  faces  the  bottom  resistance 
and  loses  its  momentum.  Thus  the  parameter  p  strongly  influences  the  jet  characteris¬ 
tics.  Reference  (39)  may  be  consulted  for  additional  details  regarding  the  jet 
solution. 

The  next  step  is  to  combine  the  jet  solution  with  the  sink  distribution 
mentioned  in  the  previous  section  to  obtain  the  jet  induced  flow  field. 

3.3  Flow  Field  Induced  By  Sink  Distribution 

Consider  a  distribution  of  sinks  representing  a  tidal  jet  along  the  positive 
x-axis  in  Figure  3.3.  In  order  to  make  the  y-axis  a  coastline  (i.e.  a  solid 
boundary)  an  image  distribution  must  also  be  introduced  as  shown. 


FIGURE  3.3  Sink  Distribution  for  a  Simple  Jet  and  Its  Image 
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I 

I 

(x,y)  is  an  arbitrary  point  in  the  flow  field,  and  the  value  of  the  stream 
function  at  (x,y)  due  to  a  sink  of  strength  q(x')  located  at  x  ■  x'  is  given  by 
(39). 


I 


'Kx.y)  - 


<Lix’) 

2ir 


tan 


-1 


y 

(x-x') 


(3.10) 


^  satisfies  Laplace's  equation  (2. 12  )  everywhere  in  the  flow  field  except  at 
x-x',  the  location  of  the  sink  itself.  If  the  sink  strength  for  the  distribution 
is  m(x)  per  unit  length  in  x  direction,  then  q(x')  -  m(x')dx'  is  the  strength 
of  an  element  dx' .  Sunning  the  stream  function  due  to  all  such  elements  of  the 
distribution,  including  the  jet  and  its  image,  the  effective  stream  function  at 
(x,y)  is 


<Kx,y) 


tan 


-1 


y 

(x-x') 


dx’ 


(3.11) 


The  preceding  superposition  holds  due  to  the  linearity  of  Eq.  (2.12).  Expressing 
the  integral  in  Eq.  (3.11)  as  the  sum  of  two  integrals. 


<Hx,y) 


m(x' ) tan 


-1 _ 

(x-x') 


dx' 


tan 


x-x 


(3.12) 


Changing  the  variable  of  integration  in  the  first  integral  on  the  RHS  from  x' 
to  -x' , 
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Mx.y)  *2 ~  J  m(-x')  tan 


x+x 


7-  dx' 


if- 


(x’)  tan  *  ^  ,  dx’ 

x-x 


(3.13) 


Due  to  symmetry  of  the  sink  distribution  about  the  y-axis. 


m(-x')  *  m(x') 


(3.14) 


and  hence  Eq.  (3.13)  becomes. 


(3.15) 


Summing  the  two  arctangent  functions  using  standard  trigonometric  identities,  the 
final  expression  for  the  stream  function  is. 


<l»(x,  y) 


(3.16) 


In  order  to  bring  the  jet  characteristics  (discussed  in  the  previous  section) 
into  the  flow  field  described  by  \l>,  the  sink  strength  m(x')  must  be  related  to 
jet  parameters.  This  is  accomplished  in  the  following  manner. 


Figure  3.4  shows  the  two-dimensional  sink  distribution  in  the  horizontal 
(x,y)  and  vertical  (x,z)  planes.  Consider  an  element  dx  of  the  distribution  located 


Volume  flow  into  the  element  demanded  by  the  sink  distribution  is 


dQ 


sink 


[m(x)dx3  h(x) 


(3.17) 


Neglecting  the  second  order  effects  of  changes  in  jet  velocity  and  width  over  a 
distance  dx,  the  volume  flow  into  dx  due  to  jet  entrainment  is  given  by 


dQent  *  (v(x,-b(x))  -  v(x,b(x))J  Jh(x)  dxj 


(3.18) 


where  v(x,-b)  is  the  component  of  jet  induced  velocity  in  the  y  direction  at  y  -  -b 
and  v(x,b)  is  similarly  defined.  Notice  that  v(x,-b)  is  positive  when  directed 
into  the  sink  whereas  v(x,b)  is  negative  under  the  same  conditions.  These  velocity 
components  at  the  edge  of  the  jet  are  referred  to  as  entrainment  velocities  in 
the  theory  of  jets.  Due  to  the  symmetry  of  the  jet  about  the  x-axis. 


v(x,-b)  -  -v (x,b) 


(3.19) 
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Since  one  needs  to  solve  the  flow  field  only  for  y  >  0  (or  y  <  0) ,  combining 
Eqs.  (3.18)  and  (3.19)  yields 

dQenC  *  -2v(x,b)h(x)dx,  y  >  0 
■  2v(x,-b)h(x)dx,  y  <  0 


(3.20) 


If  v(x,-b)  *  v  ,  so  that  v(x,b)  *  -v  ,  where  v  is  the  magnitude  of  entrainment 
e  e  e 

velocity. 


dQent  "  2v£h(x)dx  for  all  y  (3.21) 

Equating  the  entrainment  volume  flow  to  the  sink  inflow  in  Eqs.  (3.21)  and  (3.17), 

m(x)  *  2v  (x)  (3.22) 

e 

Since  h(x)  has  cancelled  out  throughout,  Eq.  (3.22)  is  applicable  for  the 
variable  bathymetry  case  also.  Now,  following  the  standard  procedure  in  turbulent 
shear  flows,  the  entrainment  velocity  is  related  to  the  centerline  velocity  of  the 
jet  through  an  entrainment  coefficient,  i.e.. 


v  (x)  =  au  (x) 
e  c 


(3.23) 


where  a  is  the  entrainment  coefficient.  Substituting  Eqs.  (3.22)  and  (3.23)  into 
Eq.  (3.16), 


<Kx 


-  ■  if 


u  (x')  tan 
c 


1  r-^_i  dx* 

l_x2_y2_x'2j 


It  is  useful  to  express  Eq.  (3.24)  in  the  nondimens ional  form 


¥(5,c) 


f  f  u(5')  tan”1!" - 2££ — ld£’ 

wJ0  U2-s2-r2J 


(3.24) 


(3.25) 


■I--'.. 
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where 


5 

V 


5 


U 


V 


JL 

b 

o 


u  b 
o  o 


(3.26) 


with  uq,  bQ  denoting 


the  inlet  throat  velocity  and  inlet  half-width  respectively. 


From  the  stream  function  in  Eq.  (3.25),  the  velocity  field  can  be  obtained 
by  the  partial  differentiation 


3j»  .:  ay 

u  -  -  ^  or  U  =  -  -r— 
3y  35 


(3.27a) 


9*  „  3¥ 

v  -  r1-  or  V  =  — 

3x  35 


(3.27b) 


giving  the  following  results: 


U(5 


.o-i-f 

,5)  -  -  f«  f 

•'n 


(g,2-g2-s2> 


(5,2-52-52)2  +  H'2c2 


d?' 


(52+C2+5’2) 


U(C’)  - 

'0  (C2-*-C2+C2)2  -  4C2C’2 


d5* 


(3.28) 


(3.29) 


Recall  from  the  discussions  in  Sec.  3.2  that  the  dimensionless  jet  centerline 
velocity  U(£)  in  Eqs.  (3.25),  (3.28),  and  (3.29)  is  given  by 
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(3.30a) 


U(C)  -  1,  C  <  Cg 


in  ZOFE  and  by 


in  ZOEF. 


(3.30b) 


The  entrainment  coefficient  a  in  the  above  equations  has  been  quoted  by 
Ozsoy  (11)  from  the  experimental  results  of  Abramovich  (13).  The  value 
of  a  is 


o  -  aj  -  0.036,  5  <  Cg  (3.31a) 

in  the  ZOFE,  and 

a  *  o2  =  0.050,  ?  >  £g  (3.31b) 

in  the  ZOEF. 


In  summary,  the  jet- induced  steady  potential  flow  in  the  presence  of  a  constant 
bathymetry  is  described  by  Eq.  (3.25)  and  Eqs.  (3.28)  through  (3.31b).  The  actual 
evaluation  and  computation  of  ¥,  U,  and  V  is  discussed  in  the  following  two  sections. 

3.4  Analytical  Treatment  Of  Solutions  For  ¥,  U,  and  V 

Eqs.  (3.25)  and  (3.28)  through  (3.31b)  show  that  they  involve  complicated 
integrals  which  cannot  be  evaluated  by  purely  analytical  means.  Prior  to  under¬ 
taking  a  fully  numerical  evaluation,  however,  it  is  useful  to  investigate  whether 
the  integrals  in  Eqs.  (3.25),  (3.28),  and  (3.29)  can  be  treated  analytically, 
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at  least  in  the  core  region,  where  the  centerline  velocity  is  a  constant  (Eq. 
(3.30a)).  This  approach  has  the  advantages  of  giving  some  insight  into  the 
solution  behavior  and  also  providing  a  check  on  the  accuracy  of  numerical  solutions. 


Using  Eqs .  (3.30)  and  (3.31),  one  can  express  the  integrals  in  Eqs.  (3.25), 
(3.28),  and  (3.29)  as  follows. 
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In  Eqs.  (3.35b),  (3.36b),  and  (3.37b),  U(0  is  given  by  Eq.  (3.30b). 

It  will  now  be  demonstrated  that  the  integrals  in  Eqs.  (3.35a),  (3.36a),  and 
(3.37a)  can  be  integrated  in  closed  form. 

First  consider  Eq.  (3.35a).  Integrating  once  by  parts  and  rearranging  the 
denominator  of  the  resulting  integrand,  it  may  be  verified  that 


(3.38a) 


with 


a  =■  1,  b  =  -2(£2-?2),  c  =  (£2+C2) 

Note  that  the  descriminant 

A  *  b2-4ac 

-  4(C2-C2)2  -  4U2+C2)2 

-  -16C2?2 

Thus,  A  <  0  for  all  £  f4  0  and  ;  t1  0 
and  A  ■  0  if  5  *  0 


(3.38b) 


(3.39) 


Notice  that  the  case  ?  *  0  is  excluded  from  the  present  considerations  because 
the  solution  is  not  expected  to  be  valid  at  the  sink  distribution  inself.  Thus, 
there  are  two  possible  cases. 
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Case  II:  £  +  0  and  £  j*  0,  i.e.,  an  interior  point  in  the  flow  field 

For  A  <  0  in  Eq.  (3.39),  defining  the  indefinite  integral  in  Eq.  (3.39a) 
as  I3, 
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It  can  be  shown  (Ref.  (40))  that. 


8ak  cos  j 


V*  -  2kC'cos  j  +  k2 

- - - - - +  - - - 

V2  +  2k£'cos  |  +  k2  4ak  sin  ~ 


[2kC *  sin  yl 

k2  -  C,2J 


which  k  =  (-)  ,  cos  e 


Combining  the  above  expression  with  Eqs.  (3.38a)  and  (3.38b),  it  can  be  shown 


simplification  that 
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(3.41) 


For  a  selected  point  5  and  c  =  1  in  the  flow  field,  it  may  be  verified 

that  Eq.  (3.41)  reduces  (after  considering  all  possible  combinations  of  signs) 


sf =  iln(4cs2 
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Now  consider  Eq.  (3.36a)  which  may  be  written  in  the  form 
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The  analytical  solutions  in  Eqs.  (3.42),  (3.45),  and  (3.47)  were  employed  to  verify 
the  correctness  and  accuracy  of  the  numerical  integration  scheme  described  in  the 
next  section. 


3.5  Development  of  Numerical  Procedure  for  Evaluation  of  S',  U,  and  V 

A  Gauss-Legendre  numerical  quadrature  procedure  (41)  was  used  to 
evaluate  the  complete  integrals  in  Eqs.  (3.25),  (3.28),  and  (3.29).  This  method 
has  an  important  advantage  over  the  straight-forward  summation  approach  in  that 
it  is  extremely  efficient  from  the  viewpoint  of  computation.  The  method  consists 
of  evaluating  the  integrand  at  selected  special  points  (the  zeros  of  Legendre 
polynomials),  multiplying  by  appropriate  weighting  factors,  and  summing  the 
resulting  values.  In  the  present  work,  sufficient  accuracy  was  achieved  by 
selecting  only  10  zeros,  whereas  a  straightforward  summation  would  require  a 
large  number  (at  least  50)  of  equi-spaced  base  points.  Therefore,  the  Gauss- 
Legendre  scheme  results  in  a  significant  saving  of  computation  time. 

The  application  of  the  proposed  numerical  scheme  requires  that  the  integrals 

in  Eqs.  (3.25),  (3.28),  and  (3.29)  be  converted  to  limits  -1  to  1.  First,  the 

upper  limit  of  08  in  these  integrals  was  replaced  by  51  5  so  that  the  integration 

s 

extends  to  50  core  lengths.  In  view  of  the  exponential  decay  of  the  jet  center- 
line  velocity  in  Eq.  (3.30b),  this  range  of  upper  limit  was  considered  sufficient. 
Next,  each  integral  was  subdivided  into  6  integrals  (0  to  Cg,  to  11  5g»  etc.,  up 
to  51  C  )  and  the  numerical  scheme  was  applied  to  each  sub-interval. 

D 
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This  subdivision  helped  to  avoid  any  loss  of  resolution  that  may  result  from 
converting  the  entire  domain  (0  to  5l£g,  or  a  large  part  thereof)  to  (-1,  1) . 
The  preceding  manipulations  are  discussed  in  detail  below. 


As  an  example  consider  the  integral  for  4*  in  Eq.  (3.25)  which  may  be 
written  as  the  sum 


*(5,0  -  £  »  (C,C) 

n  n 

n-0 


(3.48) 


where  ¥  is  given  by  Eq.  (3.35a)  and 


„  /•(10n+l)  C 

¥  U(5')  tan-1  I" - ^ - -1  d£ ' 

Q0n-9U  LC2-?2-?,2J 


(10n-9)£ 


for  n  =  1,  2,  ...  5 


(3.49) 


Similar  relations  may  be  defined  for  Eqs.  (3.28)  and  (3.29).  The  quadrature 
is  illustrated  for  Eqs.  (3.35a)  and  (3.49)  in  this  section.  First  consider 


Eq.  (3.35a)  and  define  a  new  variable 
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Changing  the  variable  of  integration  to  in  Eq.  (3.35a)  one  obtains 
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This  form  of  Eq.  (3.51)  permits  the  integration  to  be  expressed  as 
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In  Eq.  (3.52a)  0Q  denotes  a  zero  of  the  Legendre  polynomial  and  VT  is  the  cor¬ 
responding  weighting  factor,  n  is  the  total  number  of  zeros  in  the  interval 
(-1,  1)  and  its  choice  depends  on  the  accuracy  desired.  Returning  now  to  the 
remaining  integrals  in  Eq.  (3.49),  a  series  of  variables  of  integration  is 
Introduced  as  follows. 
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The  Integrals  ¥  then  become. 
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which  can  be  Integrated  to  give 
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where  the  function  is  defined  by  Eq.  (3.52b).  Following  the  above  procedures 


i 


the  integrals  in  Eqs.  (3.28)  and  (3.29)  can  also  be  evaluated  by  means  of  the  Gauss- 


Legendre  method.  Thus 
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where 
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Similarly, 
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where 
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To  summarize,  the  formulas  in  Eqs.  (3.52a),  (3.55),  (3.56a)  through  (3.56d), 

A  A 

and  (3.57a)  through  (3.57d)  were  used  for  numerical  computation  of  f,  U,  V  in 
the  flow  field.  The  actual  computational  procedure  is  described  in  the  next 
section. 


3.6  Computational  Sequence  and  Procedure 


An  efficient  computational  algorithm  was  written  to  calculate  U  and  V  at 
various  points  (£,£)  of  the  flow  field.  Typically,  the  computations  were  per¬ 
formed  in  the  following  sequence. 

(1)  For  a  given  inlet  the  friction  parameter  y  was  calculated  from  Eq.  (3.9) 

(2)  The  core  length  £g  was  computed  for  the  above  value  of  y.  Recall  that 

£  is  given  as  a  solution  of  the  condition  R(£  )  ■  0  in  Eq.  (3.1),  i.e., 
s  s 


T1  e  s  -  I2(l+aSg)  -  o 


(3.58a) 


£  was  obtained  from  this  transcendental  equation  by  a  Newton-Raphson  interaction, 
s 

For  the  special  case  of  y  ■  0,  £  is  given  by 
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(3)  As  mentioned  earlier, 10  zeros  of  Legendre  polynomials  were  chosen  and 

5'(0  )  were  computed  for  n  *  0,  1,  ....  5  from  Eqs.  (3.50b)  and  (3.53b).  These 

ni 

values  and  their  squares  were  stored. 

(4)  U{£'(0  )},  n  *  1,  2,  ...  5,  values  were  calculated  from  Eq.  (3.3) 

ni 

using  the  values  of  5'  obtained  in  the  previous  step.  These  values  were  also 
stored. 

(5)  The  computation  of  the  jet-induced  flow  field  was  then  begun  by  choosing 
a  point  (£>?)•  The  usual  practice  was  to  start  at  the  origin  and  to  proceed  along 
the  coastline  (5*0)  at  selected  intervals  of  5  up  to  5*200  (typically) .  Then,  another 
value  of  5  was  chosen  and  flow  field  computations  were  performed  at  the  same  values 

of  5.  These  calculations  were  carried  out  for  offshore  distances  of  up  to  5 1 Cs • 

The  details  of  this  procedure  are  presented  below. 

(6)  Using  the  values  of  ('(0  ),  U  (5' (9  )},  and  (5,5)  from  steps  (3), 

ni  ni 

(4),  and  (5)  respectively,  the  functions  G^,  G ^  and  G^  were  calculated  from 
Eqs.  (3.52b),  (3.56d),  and  (3.57d)  respectively,  for  each  zero  of  the  Legendre 
polynomial.  In  Eq.  (3.52b),  the  numerator  and  denominator  were  computed  separately 
for  the  argument  of  the  inverse  tangent.  The  value  of  tan  *  was  established 
by  following  an  appropriate  convention  to  account  for  its  multiplicity  w.r.t  ir. 

(7)  G^,  G2,  and  G^  were  then  multiplied  by  appropriate  constants  appearing 
in  Eqs.  (3.52a)  &  (3.55);  (3.56b)  &  (3.56c) ;  and  (3.57b)  &  (3.57c)  respectively. 

These  functions  were  summed  over  the  10  zeros  for  each  of  the  6  sub- integrals . 

Finally,  summation  was  performed  over  the  6  sub-integrals  to  obtain  the  values  of 

A  A 

U,  and  V  at  point  (5,5)  in  the  flow  field. 


The  above  numerical  scheme  was  verified  by  the  values  computed  from  Eqs. 
(3.52a),  (3.56b),  and  (3.57b)  against  the  analytical  results  of  Eqs.  (3.42), 
(3.45),  and  (3.47)  respectively.  Good  agreement  was  found. 
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4.0  POTENTIAL  FLOW  WITH  CONSTANT  BATHYMETRY  -  OTHER  CONFIGURATIONS  OF  INTEREST 

■» 

4.1  Jet  With  A  Jetty 

7- 

. 

As  outlined  in  Sec.  1.3,  the  analysis  for  a  simple  jet  presented  in  the 
previous  section  can  be  readily  extended  to  treat  other  configurations  of  tidal 
jets  and  river  mouths  as  well  (Fig.  1.2).  A  few  examples  are  given  here. 

Figure  4.1  shows  a  jet  flowing  normal  to  the  coastline  through  a  jetty  of  length 
'a'.  The  jet  is  replaced  by  an  equivalent  sink  distribution  starting  from  the 

i . 

point  x  -  a.  Also  shown  is  the  image  sink  distribution  with  respect  to  the  y-axis 
which  must  be  a  solid  boundary. 


Figure  4.1  Sink  Distribution  for  a  Jet  with  a  Jetty  and  Its  Image 

0 

The  stream  function  ifi(x,y)  at  any  point  in  the  flow  field  due  to  the  entire 
sink  distribution  can  be  obtained  following  the  method  of  Sec.  3.3.  The  final 
result  is  Eq.  (3.16)  with  the  lower  limit  ’a'  instead  of  zero.  Thus, 
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The  sink  strength  m(x)  is  also  determined  on  the  same  lines  as  in  Sec.  3.3 
(so  that  Eq.  (3.22)  holds),  except  that  the  entrainment  velocity  is  now  given 
by 


v  (x)  ■  au  (x-a) 
e  c 


(4.2) 


since  the  jet  mixes  with  the  surrounding  fluid  at  x»a.  Combining  Eqs.  (4.1), 
(4.2),  and  (3.22)  the  expression  for  is 
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Introducing  dimensionless  parameters  defined  by  Eq.  (3.26),  the  above  relation 
transforms  to 
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where  A  *  a/b 


The  corresponding  velocity  field  is  derived  from  Eqs.  (3.27a)  and  (3.27b)  to  be 
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Eqs.  (4.4),  (4.5),  and  (4.6)  completely  describe  the  flow  field  induced  by  a 
jet  with  a  jetty. 


For  convenience  of  numerical  integration,  the  expressions  for  'F ,  U,  and 
V  may  be  written  as 
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where  U(C’)  is  defined  by  Eqs.  (3.30a)  and  (3.30b)  and  a  is  given  by  Eqs.  (3.31a) 
and  (3.31b).  Numerical  evaluation  of  the  integrals  in  Eqs.  (4.7),  (4.8),  and 
(4.9)  was  performed  according  to  the  procedures  in  Secs.  3.5  and  3.6  with 
replaced  by  (C'+A)  in  Eqs.  (3.52b),  (3.56d),  and  (3.57d)  .  Computational  results 
for  this  case  are  presented  in  Sec.  5.2. 


4.2  Parallel  Jets 

To  examine  the  shelf  circulation  patterns  induced  by  multiple 
inlets,  consider  two  jets  (without  jetties)  perpendicular  to  the  coastline 

and  separated  by  a  distance  (Figure  4.2).  For  simplicity,  the  jets  considered 

are  identical  in  that  they  have  the  same  inlet  dimensions,  throat  velocities, 
and  inlet  depth.  It  is  assumed  that  the  interaction  between  the  two  jets  is 
weak,  i.e.,  the  influence  of  pressure  gradients  set  up  by  either  jet  on 
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the  other  is  neglected.  Furthermore,  the  entrainment  hypothesis  for  a  simple 
jet  must  be  modified  to  accommodate  the  case  of  two  jets.  Since  the  simple  jet 
is  symmetric  about  the  x-axis,  the  magnitude  of  the  entrainment  velocity,  v  ,  is 
the  same  on  both  jet  boundaries  (Eq.  (3.19)).  With  the  introduction  of  another 
jet,  the  fluid  in  between  the  jets  must  be  ’shared’  by  both  jets  at  their  inner 
boundaries,  whereas  an  infinite  expanse  of  fluid  is  available  on  the  outer  boun¬ 
daries  of  the  jet.  Consequently,  the  mass  entrainment  is  expected  to  be  smaller  at 
the  inner  boundaries  than  at  the  outer  boundaries.  The  jet  entrainment  was 
therefore  modified  as  follows. 


Figure  4.2  Sink  Distribution  for  Two  Parallel  Jets 


It  is  a  standard  procedure  in  the  classical  jet  theory  to  integrate 
the  depth-averaged  continuity  equation,  Eq.  (2.5),  across  the  width  of  the  jet 
to  obtain  (also  consult  Ref.  11), 


from  Eq.  (3.19).  Eq.  (4.12)  is  then  solved  simultaneously  with  the  momentum 
equation,  also  integrated  across  the  jet  width.  This  equation,  however,  does  not 
contain  any  terms  in  v^  and  need  not  be  considered  here.  If  both  sides  of  Eq. 
(4.12)  are  multiplied  by  the  density  p  (a  constant),  it  expresses  the  fact  that  the 
gradient  of  mass  flow  rate  along  the  x-axis  equals  the  total  (i.e.  both  boundaries) 
entrainment  rate  per  unit  length  across  the  x-axis.  This  is  simply  a  statement  of 
conservation  of  mass  and  it  must  hold  for  the  individual  jets  in  Figure  4.2  as 
well.  Therefore,  considering  the  lower  jet  as  an  example,  Eq.  (4.11)  is  still 
valid,  but  the  entrainment  velocities  v(x,-b)  and  v(x,b)  are  different.  At  this 
point,  due  to  the  lack  of  experimental  entrainment  data  for  the  two  jets  configura¬ 
tion,  it  is  assumed  that  the  total  mass  flow  rate  is  the  same  as  that  for  a  simple 
jet.  In  other  words,  the  drop  in  entrainment  on  the  inner  boundary  equals  the 
increase  in  entrainment  on  the  other  boundary.  Thus,  in  Eq.  (4.11), 

v(x,-b)  »  v  +6  (4.13a) 

e 

v(x,b)  =  -(v  -6)  (4.13b) 

e 

where  vg  is  the  entrainment  velocity  for  a  simple  jet.  Substituting  Eqs.  (4.13a) 
and  (4.13b)  into  Eq.  (4.11),  Eq.  (4.12)  is  recovered.  The  change  in  entrainment 
velocity  6  obviously  depends  upon  the  relative  strengths  (such  as  inlet  velocity, 
width,  etc.)  of  the  two  jets.  In  the  example  of  Figure  4.2,  however,  the  jets  are 
assumed  to  have  equal  strengths  for  simplicity. 
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Consider  the  two  jets  and  their  images  as  shown  in  Figure  4.2.  For  the 
lower  jet,  denoted  by  subscript  1,  Eq.  (3.25)  gives  the  expression  for  the  stream 
function. 


*i«,n) 


U(E')  tan 


(4.14) 


A  similar  expression  can  be  written  for  the  upper  jet  (i.e.,  jet  2),  where  5  must  be 
replaced  by  (£-L)  since  the  x-axis  is  located  along  the  centerline  of  jet  1.  Thus, 


*2(5,n) 


U(C')  tan 


-1 


25U-L) 


C2-(5-L)2-E’2 


(4.15) 


Notice  that  both  Eqs.  (4.14)  and  (4.15,/  contain  the  same  centerline  velocity 
U(5)  since  the  two  jets  are  identical.  The  combined  effect  of  two  jets  at  any 
point  (£,C)  in  the  flow  field  is  obtained  by  the  superposition, 

’  -  *1  +  *2  (4.16) 


or 


(4.17) 


The  jet- induced  velocities,  U  and  V  may  also  be  obtained  by  similar  super¬ 
position  or  by  carrying  out  the  differentiation  in  Eqs.  (3.27a)  and  (3.27b). 
The  result  is 
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Eqs.  (A.  17),  (4.18),  and  (4.19)  describe  the  flow  field  induced  by  two  parallel 


)ets,  subject  to  the  various  assumptions  outlined  in  the  beginning  of  this  section 


It  Is  not  necessary  to  further  simplify  the  above  equations,  since  each  integral 


can  be  evaluated  separately  using  the  numerical  procedures  in  Secs.  3.5  and  3.b 


The  results  are  given  in  Sec.  5.3 


4.3  Jet  Non- Per pend  leu la r  to  Coastline 


The  simple  Jet  in  Fig.  4.1  shown  normal  to  the  coastline  probably  represents 


an  idealization.  In  reality,  the  tidal  Jet  Is  more  likely  to  be  depleted  as 


As  before,  the  jet  Is  replaced  by  a  sink  distribution  and  by  an  image  distribution 
in  order  to  make  the  coastline  a  solid  boundary.  The  coordinate  systems  x,  y 
(for  real  jet)  and  x,  y  (for  image  jet)  have  been  introduced  for  convenience  of 
mathematical  derivation.  Once  again,  the  entrainment  velocity  induced  by  the  jet 
must  be  modified  as  in  the  previous  section.  In  the  present  case,  the  entrainment 
at  the  jet  boundary  nearer  to  the  coastline  is  expected  to  be  less  than  the 
entrainment  on  the  other  boundary.  The  same  principle,  as  outlined  in  Sec.  4.2, 
is  extended  to  the  configuration  in  Fig.  4.3,  and  the  total  mass  entrainment  is 
assumed  to  equal  that  for  a  simple  jet.  Therefore,  Eq.  (4.12)  and  the  Jet  solution 
given  in  Sec.  3.2arestill  valid.  The  jet  induced  flow  field  is  determined  as 
follows. 


The  dimensional  stream  function  ip  due  to  the  real  jet  (w.r.t.  x,  y  co-ordinates) 
is  given  by  an  equation  similar  to  Eq.  (3.11). 


*r(x,y)  -  ± 


/ 


m(x')  tan 


-1 


(x-x’) 


dx' 


(4.20) 


The  stream  function  due  to  the  image  Jet  is 


*im  <*•*>  ■  27 


/ 


m(x')  tan 


-1 


(x— x 1  )j 


bx' 


(4.21) 


Eqs.  (4.20)  and  (4.21)  must  be  transformed  to  x,  y  coordinates  using  the  relations 


x  ■  x  cos  $  +  y  sin  <J> 


(4.22a) 


y  ■  -x  sin  •C  +  y  cos  $ 


and 
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x  *  x  cos  P  -  y  sin 


(4.23a) 


y  *  x  sin  <J>  +  y  cos 


(4.23b) 


Substituting  Eqs.  (4.22)  and  (4.23)  into  Eqs.  (4.20)  and  (4.21)  respectively, 

and  superposing  i|>  and  i|i  ,  the  stream  function  in  x,y  coordinates  is 
lm  r 


'Kx.y)  » 


2tt 


y  co_s.t_...  d°, 


sin  p  -  x* 


2tt 


f  .<;■>  tan'4  -x  sl"  i  ■+ 

J  lx  cos  p  -  y 

■>00  I 

f  ,(?)  tan'1  [  --  -x-31-"-  ±  1  ■9°S-J’.--]d;> 

g*'  lx  cos  p  +  y  sin  $  -  x'J 


(4.24) 


Changing  the  variable  of  integration  from  x'  to  -x'  in  the  first  integral  of  the 
above  equation  and  making  use  of  the  symmetry  of  the  sink  distribution  w.r.t.  the 
y-axis,  i.e.,  m(x’)  »  -  m(-x’),  the  result  is. 


P(x,y) 


I  <m(x’)<tan' 


tan 


dx'  (4.25) 


The  arctangent  functions  in  Eq.  (4.25)  are  summed  using  standard  trigonometric 
formulas  and  the*  sink  strength  m  is  substituted  from  Eqs.  (3.22)  and  (3.23)  to  obtain 


>Kx,y)  *  ~ 


/« 

*>’)  tan-1f- 

u  l> 


2xy  -  2xx'  sin 


u  (x’)  tan 

lx2-y2-?2  +  2x'y  sin  p 


J  dx' 


(4.26) 


Notice  that  Eq.  (4.26)  reduces  to  Eq.  (3.24)  for  a  simple  jet  when  p  ■  0. 
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Replacing  the  dummy  variable  of  integration  x'  by  x'  without  loss  of  generality, 
and  non-dimensionalizing  using  Eqs.  (3.26),  the  above  equation  becomes, 


Y(€,0  -  J 


M  OO 

I  U(C')  tan’1  [  26.g  -  sin  f  1 

/  Ls2-C2-C* 2  +  25’c  sin  4>J 


dC' 


(4.27) 


The  velocity  field  is  obtained  by  differentiating  Eq.  (4.27)  according  to  Eqs. 
(3.27)  to  get 


U(£,C) 


V(5,C) 


NjU.C,?',*) 

D  U,C,5',*)  dC 

N2(C,C,?',<fr) 

D  (C, C,S\4>)  dC' 


(4.28) 


(4.29) 


where , 

N1(S,?,f-\<j>)  -  (C,2-C2-C2)  +  2C'  sin  4)  (?-£'  sin  <(>)  (4.30a) 

-  (C-C  sin  <J>)  (52+?2+5,2-25'5  sin  4.)  (4.30b) 

D  (5,C,£\$)  -  452(c-5'  sin  4>)2  +  (52-?2-€'2  +  2C'C  sin  4>)2  (4.30c) 

It  may  be  verified  that  Eqs.  (4.28)  and  (4.29)  together  with  Eqs.  (4.30),  revert 
to  Eqs.  (3.28)  and  (3.29)  when  4>  =  0. 

The  flow  field  induced  by  a  jet  which  is  not  perpendicular  to  the  coastline 
is  completely  described  by  Eqs.  (4.27),  (4.28),  and  (4.29)  supplemented  by  Eqs. 
(4.30).  As  in  Secs.  4.1  and  4.2,  numerical  results  were  obtained  by  the  procedures 
discussed  in  Secs.  3.5  and  3.6,  and  the  results  are  presented  in  Sec.  5.3. 


5.0  CONTINENTAL  SHELF  CIRCULATION  RESULTS 


5. 1  The  Simple  Jet 

The  flow  field  induced  by  a  tidal  jet  modeled  as  a  simple  jet  is  described 
by  Eqs.  (3.25)  and  (3.28)  through  (3.31).  Numerical  evaluation  of  V,  U,  and  V 
in  these  equations  was  carried  out  following  the  procedures  of  Secs.  3.5  and 
3.6.  The  computational  results  are  presented  in  Figs.  5.1  through  5.5. 

Figure  5.1  shows  the  streamline  pattern  for  the  ideal  case  of  a  jet  with 
zero  bottom  friction  (p*0) .  Streamlines  have  been  plotted  only  for  C  >  0  since 
the  flow  field  is  symmetric  about  the  5*0  axis.  The  lines  of  constant  V 
were  obtained  by  inverse  linear  interpolation  through  the  values  calculated  by 
Eq.  (3.25).  Due  to  the  relatively  slow  linear  growth  of  the  bottom-frictionless 
jet,  as  opposed  to  the  exponential  growth  with  bottom-friction,  the  jet  induces 
circulation  over  a  large  portion  of  the  continental  shelf  (up  to  400  inlet  half¬ 
widths  offshore  in  Fig.  5.1).  The  flow  direction  is  along  the  streamlines  into 
the  jet  with  a  small  onshore  component  which  vanishes  along  the  coastline. 

The  case  of  a  bottom  frictional  jet  is  shown  in  Fig.  5.2  for  p  *  0.05.  This 
value  is  expected  to  be  representative  for  a  typical  inlet  (such  as  Redfish 

Pass,  Florida  with  b  /h  =  15)  with  a  Darcy-Weisbach  friction  factor  f  *  0.02  to  0.08. 
o  o 

A  comparison  of  Figures  5.1  and  5.2  reveals  that  bottom  friction  in  the  jet  has 
a  significant  effect  on  the  shelf  circulation  pattern.  The  fluid  enters  the 
jet  almost  radially  as  shown  by  the  streamlines  and  the  influence  of  the  jet 
is  confined  to  a  small  distance  offshore.  This  is,  of  course,  due  to  the 
rapid  exponential  growth  of  the  jet  near  the  inlet  as  shown  in  Figure  3.2. 
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Figure  5.3  displays  the  jet  induced  flow  pattern  for  a  very  large  value  of  the 
friction  parameter  (u  »  0.10).  Convergence  of  entrainment  streamlines  toward 
the  inlet  is  very  pronounced  in  this  case.  Since  all  of  the  streamline  patterns 
discussed  were  generated  for  jets  replaced  by  sink  distributions  along  their 
axes,  the  solutions  shown  are  not  applicable  within  the  physical  jet  boundaries 
which  have  been  superimposed  in  Fig.  5.1,  5.2,  and  5.3. 

An  important  aspect  of  the  jet-induced  flow  fields  described  so  far  is  the 
variation  of  the  alongshore  current.  Since  entrainment  due  to  the  tidal 
jet  is  a  persistent  phenomenon  during  every  ebb  tide,  the  alongshore  current 
can  be  a  dominant  cause  of  sediment  transport  toward  the  inlet  from  the  adjacent 
beaches.  Figure  5.4  shows  the  variation  of  dimensionless  alongshore  velocity 
with  distance  from  the  inlet.  The  upper  three  curves  correspond  to  p  *  0,  0.05, 
and  0.1  for  a  simple  jet  (l.e.,  without  jetty,  A  ■  0)  and  the  bottom  curve 
corresponds  to  a  jet  with  a  jetty.  The  latter  is  included  for  the  purpose  of 
comparlsion  and  will  be  discussed  in  the  next  section.  It  is  seen  that  near 
the  inlet  the  alongshore  velocity  is  roughly  3.6%  of  the  inlet  throat  velocity 
for  all  friction  parameters.  This  particular  value  is  a  direct  result  of  the 
entrainment  coefficient  oj  ■  0.036  in  the  core  region  of  the  jet  (Eq.  (3.31a)). 

As  one  moves  away  from  the  inlet,  along  the  coastline,  the  alongshore  velocity 
for  Che  bottom-frictionless  jet  decreases.  With  the  introduction  of  friction 
(p  ■  0.05)  the  reduction  in  velocity  with  distance  from  the  inlet  is  quite 
rapid  as  shown  by  the  top  two  curves  in  Figure  5.4.  Further  increases  in  friction 
accentuates  this  drop-off  of  the  alongshore  current. 

It  is  also  interesting  to  consider  the  behavior  of  the  u-component  of  velocity 

A 

in  the  flow  field.  Fig.  5.5  shows  the  variation  of  U  with  offshore  distance  for 


y.fkyiit.' 
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different  locations  of  ?  from  the  inlet.  Since  the  u-  velocity  is  zero  at  £  =  0 
and  since  it  must  also  be  zero  far  away  from  the  jet,  the  curves  of  U  exhibit 
a  peak  as  seen  in  the  figure.  The  magnitude  of  U  at  the  peak  decreases  with  the 
increase  in  the  distance  £  from  the  jet.  This  is  expected  because  fluid  particles 
closer  to  the  jet  experience  a  greater  influence  due  to  the  jet.  It  is  also 
seen  tnat  the  point  at  which  the  maximum  of  U  occurs  moves  further  offshore  as 
cr.e  moves  away  from  the  jet. 

5.2  Jet  With  A  Jetty 

Streamlines  for  this  case  were  obtained  through  numerical  evaluation  of 
Eqs.  (4.4),  (4.5),  and  (4.6).  The  results  are  shown  in  Figs.  5.6,  5.7,  and  5.8 
for  p  **  0.1  and  dimensionless  jetty  lengths  A=  2,  5,  and  10  respectively.  These 
figures  should  be  compared  with  the  earlier  Fig.  5.3  which  displays  streamlines 
for  a  jet  without  jetty  and  p  *  0.1.  While  Figs.  5.6  through  5.8  display  the 
radially  converging  streamlines  associated  with  the  high  value  of  friction  para¬ 
meter  (p  *  0.1),  a  significant  curvature  of  streamlines  close  to  the  jetty  is 
observed.  In  other  words,  fluid  particles  adjacent  to  the  jetty  are  entrained 
into  the  jet  along  streamlines  curving  toward  the  discharge  point  of  the  jetty. 

The  variation  of  alongshore  current  is  especially  important  in  this  case 
since  accretion  of  material  near  the  mouths  of  jetties  is  a  common  observation.  Figure 
5.9  contains  a  plot  of  dimensionless  alongshore  curren£.against  dimensionless  off¬ 
shore  distance  for  different  values  of  jetty  length  and  p  =  0.01.  The  curve  for 
a  jet  without  a  jetty  has  been  reproduced  from  Figure  5.4  and  is  included  for 

A 

comparison.  First,  as  expected,  the  alongshore  velocity  V  goes  to  zero  at 
the  jetty  (4  •  0)  since  it  is  normal  to  a  solid  boundary.  At  some  distance 

A 

from  the  jetty,  V  exhibits  a  maximum  and  then  decays  to  zero.  The  peak  value 
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of  alongshore  current  drops  sharply  with  increasing  jetty  length.  In  fact, 
introduction  of  a  jetty,  which  is  only  two  inlet  half -widths  long  (A=2) ,  reduces 
the  maximum  V  by  a  factor  of  1/2.  Moreover,  the  distance  from  the  jetty  at 
which  the  maximum  velocity  occurs  increases  with  each  increase  in  jetty  length. 

5.3  Parallel  Jets 

The  streamline  pattern  for  two  parallel  jets,  obtained  from  Eq.  (4.17)  is 
shown  in  Figure  5.10.  Both  jets  were  considered  bottom-frictionless  since  this 
allows  a  slow  linear  growth  (as  opposed  to  an  exponential  growth  for  a  bottom- 
frictional  jet)  before  the  boundaries  of  the  jets  intersect.  The  jets  in  Figure  5.10 
are  separated  by  a  distance  of  50  inlet  half-widths,  which  is  typical  for  many 
adjacent  inlets.  The  flow  field  is  shown  for  y  >  0  only  since  it  is  expected  to 
be  symmetric  about  the  centerline  between  the  two  jets  (i.e.  c  *  25).  The 
streamlines  between  the  jets  are  parallel  to  the  centerline  away  from  the  inlets 
and  display  strong  curvature  due  to  the  entrainment  effect  on  either  side  of 
the  centerline.  The  streamlines  in  the  outer  region  (e  >  50  and  5  <  0  by  symmetry) 
resemble  closely  the  pattern  for  a  simple  jet  (Fig.  5.1)  It  is  seen  that 
the  coastline  and  the  centerline  between  the  jets  correspond  to  a  single  stream¬ 
line  because  the  velocity  component  normal  to  these  lines  is  zero.  Thus,  the 
centerline  (£  ■  25)  acts  as  a  dividing  streamline.  In  the  zone  U  <  ?  <  25,  the 
fluid  is  entrained  into  the  lower  jet  (labeled  Jet  1)  and  in  the  zone  25  <  C  <  50, 
the  fluid  is  entrained  into  the  upper  jet  (labeled  Jet  2). 

The  foregoing  results  obviously  correspond  to  the  mathematical  idealization 
of  jets  acting  as  distributed  sinks  (Fig.  4.2).  In  reality,  the  physical  dimen¬ 
sions  of  the  jets  must  be  considered.  In  Fig.  5.10,  the  jet  boundaries  are 
superimposed  over  the  streamline  pattern  to  delineate  the  region  of  applicability 
of  the  theoretical  solution.  For  ?  >  50,  this  region  contains  the  flow  field 
between  the  jet  boundary  and  the  coastline,  and  for  0  <  5  <  50,  the  region 


4, 0'VkVH,'^ 


includes  the  fluid  bounded  by  the  two  jet  boundaries  and  the  coastline. 

The  variation  of  alongshore  current  in  the  present  case  is  most  interesting 
(Fig.  5.11).  For0<£<25,  the  vertical  component  velocity  V  is  directed  into  the  lower 
jet  (Fig.  5.10)  and  hence  negative.  For  25  <  C  <  50,  V  is  directed  into  the 
upper  jet  and  hence  positive.  At  the  centerline  between  the  two  jets,  i.e., 

A  A 

?  »  25,  V  equals  zero.  This  variation  in  V  is  antisymmetric  since  the  two  jets 
are  identical.  V  changes  from  -0.01  near  5  *  0  to  +0.01  near  C  =  50,  passing 
through  zero  at  c  =  25.  Notice  that  the  magnitude  of  V  near  the  jets  is  much 
lower  than  that  for  a  single  jet  (0.01  compared  to  0.036  for  the  latter).  This 
is  due  to  the  'sharing'  of  fluid  between  the  two  jets  as  pointed  out  in  Sec. 

4.2.  When  one  crosses  the  jet  at  £  =  50  moving  along  the  coastline,  the  along¬ 
shore  current  jumps  discontinuously  to  V  »  -  0.062  which  is  much  greater  than  the 
value  0.036  for  a  simple  jet.  It  may  be  recalled  from  Sec.  4.2  that,  for  the 
parallel  jet  situation,  the  entrainment  on  the  inner  boundaries  decreases 
whereas  the  entrainment  on  the  outer  boundaries  increases  such  that  the 
total  mass  (or  volume)  entrained  is  the  same  as  that  for  a  single  jet.  This 
inequality  in  entrainment  on  the  boundaries  of  each  jet  is  a  consequence  of  mass 
conservation.  The  discontinuity  in  V  across  the  jet  at  £  =  50  (since  a  solution 

A 

does  not  exist  at  the  jet  axis  itself)  from  |v|  =  +0.01  to  0.062  preserves  the 
total  mass  entrainment  of  0.072  for  a  simple  jet. 

The  above  results  indicate  that,  for  two  adjacent  jets,  the  alongshore 
current  in  the  outer  region  is  much  stronger  than  the  current  in  between  the  jets. 
Therefore,  there  should  be  an  increased  transport  of  suspended  sediment  toward 
the  two  inlets  along  the  beaches  external  to  the  two  jets.  On  the  other  hand, 
transport  of  material  toward  the  inlets  in  between  the  two  jets  should  be  com¬ 
paratively  much  smaller. 
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5.4  Jet  Non-Perpendicular  to  Coastline 

Eqs.  (4.27),  (4.28),  and  (4.29)  were  evaluated  numerically  to  obtain  the 
stream  function  ¥  and  the  velocities  U  and  V  in  the  flow  field.  Fig.  5.12 
shows  the  streamline  pattern  for  a  simple  jet  (with  y  *  0.05)  which  has  been 
rotated  anticlockwise  through  15°.  The  flow  field  is  quite  similar  to  that 
in  Fig.  5.2  for  a  simple  jet  perpendicular  to  the  coastline.  The  asymmetrical 
nature  of  the  jet  induced  flow  is  obviously  due  to  the  non-perpendicularity  of 
the  jet. 


The  alongshore  current  is  plotted  in  Fig.  5.13.  As  an  observer  moves  along 
the  coastline,  starting  at  C  *  -»  (i.e.  away  from  the  jet  on  the  negative  axis  in 
Fig.  4.3),  the  velocity  V  Increases  gradually  from  0  to  0.05  for  £  just  less 
than  zero.  When  the  observer  crosses  the  ?  ■  0  axis  on  the  other  side  of  the 
jet  (Fig.  4.3),  the  V  changes  discontinuously  to  -0.022.  Thus,  the  entrainment 
is  greater  on  the  lower  boundary  of  the  jet  and  smaller  on  the  upper  boundary. 
This  is  in  accordance  with  the  mass  conservation  considerations  discussed  in 
Sec.  4.3.  Therefore,  the  alongshore  current  on  the  obtuse  angle  side  of  a 
non-perpendicular  jet  is  stronger  than  the  current  on  the  acute  angle  side. 
Consequently,  there  should  be  increased  transport  of  sediment  toward  the  inlet 
on  the  obtuse  angle  side  compared  to  the  acute  angle  side. 
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DISTANCE  ALONGSHORE 


6.0  VARIABLE  DEPTH  AND  UNSTEADY  EFFECTS  -  PRELIMINARY  CONSIDERATIONS 

6. 1  Potential  Flow  with  Variable  Depth 

The  work  presented  in  Sections  3.0,  4.0,  and  5.0  was  concerned  with  jet- 
induced  circulation  treated  as  a  steady,  potential  flow  over  a  constant 
bottom.  In  this  section,  preliminary  formulations  are  derived  to  include  the 
effects  of  variable  bathymetry,  unsteadiness,  and  the  Earth's  rotation.  The 
variable  depth  problem  is  considered  first.  In  Sec.  2.2  it  was  shown  that 
steady  potential  flow  over  arbitrary  bottom  topography  is  governed  by  Eq.  (2.8) 
in  terms  of  a  velocity  potential  and  by  Eq.  (2.11)  in  terms  of  a  stream  function. 
The  latter  formulation  was  employed  in  the  constant  bathymetry  case  which  is 
governed  by  Laplace's  Eq.  (2.12).  General  solutions  were  obtained  by  replacing 
the  tidal  jet  by  an  equivalent  sink  distribution  and  superimposing  the  elementary 
solution  (Eq.  (3.10))  for  each  sink.  Had  this  problem  been  worked  in  terms  of  a 
velocity  potential  instead  of  the  stream  function,  the  corresponding  elementary 
solution  due  to  a  sink  located  at  the  origin  (x'=y'=0)  would  be 

4>(r)  -  ^  £n  r  (6.1) 

where  r=^x^+y^.  The  potential  function  in  Eq.  (6.1)  has  the  important  pro¬ 
perty  that  it  posesses  a  logarithmic  singularity  at  the  origin.  This  is  reflected 
in  the  form  of  an  indeterminate  argument  and  resulting  multiplicity  of  the 
arctangent  function  at  the  origin  in  the  stream  function  formulation  (Eq.  (3.10)). 
These  behaviours  of  the  potential  and  stream  functions  are  due  to  the  fact 
that  Laplace’s  Eqs.  (2.9)  and  (2.12)  are  not  satisifed  at  the  sink  itself. 


Consider  now  the  governing  equations  for  the  variable  depth  problem, 
viz.,  Eqs.  (2.8)  and  (2.11).  As  discussed  in  Ref.  (42),  it  has  been  demonstra¬ 
ted  by  Hadamard  (43)  that  for  the  general  elliptic  equation  in  the  plane, 
such  as  Eq.  (2.8),  an  elementary  solution  containing  a  singularity 
exists  and  it  is  unique.  A  method  for  solving  Eq.  (2.8)  analytically  is  also 
given  in  Ref.  (42).  Therefore,  the  velocity  potential  approach  was  selected  over 
the  stream  function  approach  to  analyze  the  variable  depth  problem.  The  next 
three  sections  deal  with  the  derivation  and  application  of  solutions  of  Eq.  (2.8) 


6.2  General  Elementary  Solution  for  Arbitrary  Depth 


Eq.  (2.8),  repeated  below  for  convenience,  is 


h  h 

V2<(>  +  +  +  -£■  <t> 

h  x  h  y 


(6.2) 


If  Eq.  (6.2)  is  non-dimensionalized  using  Eqs.  (3.26),  the  result  is 


H  H 

V24>  +  — i  <j,  +  — 2.  $ 

H  5  H  £ 


(6.3) 


which  is  identical  in  form  to  Eq.  (6.2).  Thus,  it  is  sufficient  to  consider 
Eq.  (6.2)  and  treat  the  quantities  $,  h,  x,  and  y  as  dimensionless.  Eq.  (6.2) 
may  be  written  in  the  general  form. 


L(u)  *  0 


(6.4a) 


where  the  linear  operator  L  is  given  by 


L(u)  =  V2u  +  2au  +  2bu  +  cu 
x  y 


(6.4b) 
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First,  an  elementary  solution  of  the  generalized  equation  (6.4)  is  sought  in  the 
form 


u 


n=0 


U  D 
n 


n+v 


(6.5) 


where 


D  =  r2 

U  -  U  (r ,  6) 
n  n  ’  ' 

with 

r  *  V  x2+y2 

and 


(6.6a) 

(6.6b) 


(6.6c) 


The  series  in  Eq.  (6.5)  is  equivalent  to  the  Method  0f  Frobenius  for  ordinary 
differential  equations.  The  constant  v  is  to  be  suitably  chosen  and  the  co¬ 
efficients  Un  are  analytic. 


The  details  of  the  substitution  of  Eq.  (6.5)  into  Eq.  (6.4)  and  determination 
of  the  solutions  are  given  in  Ref.  (42).  These  derivations  are  tedious  and  only  the 
final  results  are  presented  here.  It  can  be  proved  that  there  are  two  linearly 
independent  solutions  of  Eq.  (6.4).  The  first  solution  is 


w 


U  -  U  £  -J  Dn 
1  °  n*0  n! 


(6.7a) 


where. 
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w  =  1  and 
o 


w  » - —  (*  L,  (w  ,)rn  1 

"  4i"  J  1 


(6.7b) 


T  ,  \  32w  ,  32w  .  _  3w  „  3w 

L,  (w)  =  -  -t-  -  +  2a,  —  +  2b,  —  +  c,w 

1  3x2  3y2  1  3x  1  3?  1 


(6.8a) 


,  3U  ,  3U  ,  3U 

_  .  1  O  ,  ,  ,  1  O  lo 

al  +  U  3x  ’  bl  =  b  +  U  3y  ’  cl  =  U  3y 
o  o  o 


(6.8b) 


In  Eqs.  (6.7)  and  (6.8),  Uq  is  given  by 


U  =  A  exp 
o 


{-/•■ 


(6.9a) 


where  the  constant  A  may  depend  on  9  and 


a  =  ax  +  by 


(6.9b) 


A  second  solution  of  Eq.  (6.4)  is  given  by 


•  f  a  (F(vfl)  W 

uo  £ jrc^firj 


a  Dn 


(6. 10a) 


where 


W  *  1 ,  W  - - f  L,  (W  )  rn+v_1 

o  n  4rn+v  J  1  n-1 


(6.10b) 


In  Eq.  (6.10a),  T  denotes  the  Gamma  function.  It  is  noticed  that  the  first 
solution  in  Eq.  (6.7a)  is  analytic  in  a  neighborhood  of  the  origin.  The  second 
solution  in  Eq.  (6.10a)  has  a  term  containing  a  logarithmic  singularity  and  a 


86 


AD-A075  158  TETRA  TECH  INC  JACKSONVILLE  FL  F/G  8/3 

CONTINENTAL  SHELF  CIRCULATION  INDUCED  BY  TIDAL  JETS  -  AN  ANALYT— ETC 
OCT  79  R  B  TAYLOR  #  P  B  JOSHIO  N00014-78-C-0693 


UNCLASSIFIED  TETRAT-3249  NL 


2  OF  2 

^076158 

■ 

b 

■ 

|§ 

H 

i§ 

i 

END 

DATE 

FILMED 

11-79 

DOC 

term  which  is  analytic  in  a  neighborhood  of  the  origin.  This  solution 
is  the  required  elementary  solution  for  the  general  equation  (6.4).  For 
details  regarding  the  convergence  of  the  above  solutions  Res.  (42)  and  (43) 

should  be  consulted. 

In  the  next  section,  the  general  solutions  discussed  heretofore  are 
applied  to  the  special  case  of  an  exponentially  sloping  bottom. 


6.3  Derivation  of  Elementary  Solution  for  an  Exponentially  Sloping  Bottom 


Consider  a  special  case  of  Fig.  L.3  in  which  the  ocean  floor  slopes  exponen¬ 
tially  from  the  coast  in  the  x-direction,  i.e.. 


h(x)  *  e 


2Bx 


(6.11) 


where  h  is  dimensionless  depth  which  equals  unity  at  the  inlet  (x“0)  and  8  is 
a  non-dimensional  slope  parameter  .  In  fact,  26  equals  the  slope  at  the  origin,  since 

28x 


h  -  28e 
x 


(6.12a) 


and 


h  •  0 

y 


Therefore,  in  Eqs.  (6.4),  in  conjunction  with  Eq.  (6.2), 


2*  '  T  '  26 


or 


(6.13a) 


a  ■  6  ■  constant 


(6.14a) 
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(6.14b) 


I 

I 

t 

[ 

r. 


and 


.  n 

»-T-f  50 


Substituting  Eqs.  (6.14)  into  Eqs.  (6.9)  and  using  Eqs.  (6.6),  it  can  be  shown 
that 


U  -  Ae"*X 
o 


(6.15) 


Now,  from  Eqs.  (6.7b)  and  (6.8), the  coefficients  w  can  be  evaluated  as 

n 

,  62  84  86 
"o'1'  "l  ■  T'  "2  ■  32  •  vm  •  ••• 


or 


8 


2  n 


22n  n! 


(6.16) 


Introducing  Eq.  (6.16)  into  Eq.  (6.7a),  the  first  solution  is 

»2«2 /Ain 


.  Ae-8x  £  I8f r^ 
n*0  (n ! ) 2 


(6.17) 


In  the  above  equation,  A  is  a  pure  constant  Independent  of  6.  In  fact,  it  can 
be  proved  that  A  must  be  a  constant  in  order  to  preserve  the  analytlclty  of  u^ 
at  the  origin  and  to  preserve  symmetry  of  the  resulting  velociy  field  about 
0*0. 


To  obtain  the  second  solution  in  Eq.  (6.10a),  the  coefficients  are 


evaluated  from  Eq.  (6.10b).  The  result  is 


"  lf  ^  “  4(\H-1)  *  W2  "  16(v+l)  (v+2)  » 


g6 

and  ^  "  64(v+l)  (v+2)  (v+3)  ' 
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or 


S2n 


22?  FI  (v+n-r) 
r-0 


e20  r(v+l) 

2 211  r  (v+n+l) 


(6.18) 


Substituting  the  series  (6.18)  into  Eq.  (6.10a),  the  required  elementary  solution 
is. 


u 


2 


u.  Zn  D  +  U 
1  o 


09 

£ 

n»0 


g2n 

22n 


3_ 

3v 


Ir(u-H) 

|  r (v+n+l) 


v«0 


(6.19) 


When  the  differentiation  in  the  above  equation  is  carried  out,  and  u^  and  are 
substituted  from  Eqs.  (6.17)  and  (6.15)  respectively,  the  result  is 


(82r2/4)n  tn  r  _  ^  l  (g2r2/4)2m' 

(n!)2  m^f  p*l  p  (m!)2 


(6.20) 


u2  is  the  elementary  singularity  solution  for  a  sloping  bottom.  It  can  be 
shown  that  the  series  in  Eqs.  (6.17)  and  (6.20)  are  absolutely  and  uniformly 
convergent.  The  solutions  given  by  Eqs.  (6.17)  and  (6.20)  appear  complicated, 
but  fortunately,  the  summations  for  u^  and  u2  tend  in  the  limit  to  the  well- 
known  Bessel  functions.  It  may  be  verified  from  Ref.  (44)  that 

uL  -  A  e*Sx  lo(8r)  (6.21) 

u2  -  -2Ae'0x  [xo(8r)  +  ygr)  In  (8eE/2)j  (6.22) 

where  1q  and  are  zeroth  order  modified  Bessel  functions  of  the  first  kind 
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and  second  kind  respectively.  E  is  Euler's  constant  equal  to  0.577215664... 
It  can  be  demonstrated  by  direct  substitution  that  Eqs.  (6.21)  and  (6.22) 
satisfy  the  governing  Eq.  (6.2). 

In  Eq.  (6.22),  K  (8r)  behaves  like  ta  r  as  r-*-0 ,  and  hence  u.  is  the 
o  t 

required  elementary  solution  for  $  in  Eq.  (6.2).  Thus,  for  an  exponentially 
sloping  bottom,  the  velocity  potential  in  polar  co-ordinates  is 


*(r,8)  -  -2Ae'8rCOs6  [^(Br)  +  IQ(Br)  £n  (BeE/2)j 


Corresponding  velocity  components  are  given  by  the  derivatives. 


a  .21 

r  3r 


.  .1  21 

8  r  38 


(6.23) 


(6.24a) 

(6.24b) 


Differentiating  Eq.  (6.23)  one  obtains 

u  (r,8)«2A Be  ®rcos®  f{K.(6r)  +  K  (Br)  cos8} 
r  l  l  o 

-  txi  (BeE/2)  {I^Br)  -  Io(8r)  cos8}]  (6.25) 

uQ(r,8)  «  -2A8  sine  e_6r  Cos6  [^(Br)  +  £n(BeE/2)lo(Br)j  (6.26) 


Note  that  u  is  an  even  function  of  8,  i.e.,  u  (r,9)  »  u  (r,-9),  and  uQ  is  an  odd 
function,  i.e.,  UQ<r,9)  »  -u^(r,-8).  In  Eq.  (6.25),  1^  and  are  first  order  modi¬ 
fied  Bessel  functions  of  the  first  and  second  kind  respectively.  It  should  be  empha¬ 
sized  that  Eqs.  (6.23),  (6.25),  and  (6.26)  describe  the  potential  flow  field  about  a 
point  singularity  for  exponentially  varying  depth.  To  obtain  solutions  for 
outer  shelf  circulation  one  must  relate  the  jet  solution  for  exponential  bathy¬ 
metry  to  a  distribution  of  singularities  represented  by  Eq.  (6.23),  (6.25),  and 

(6.26).  This  precedure  is  outlined  in  Sec.  6.4 
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Finally,  Che  constant  A  in  Eqs.  (6.23),  (6.25),  and  (6.26)  must  be  given 
a  physical  interpretation  by  identifying  it  with  the  strength  of  the  singularity. 
By  analogy  with  the  elementary  solutions  for  the  constant  bottom  case  (Sec.  3.1), 
the  singularity  may  be  referred  to  as  a  'modified'  sink.  Let  'q'  be  the 
strength  of  the  modified  sink  at  r»0.  Thus,  q  equals  the  volume  of  fluid 
vanishing  per  unit  time  at  the  sink  per  unit  length  (measured  normal  to  the  plane 
of  motion).  Then,  q  also  equals  the  volume  of  fluid  flowing  into  any  curve 
enclosing  the  sink.  For  simplicity,  let  the  curve  be  a  circle  of  vanishingly 
small  radius  e  with  center  at  r*0.  Therefore, 

q  -  (f)  V.  er  dX  (6.27) 

circle 

e 


where  dX  is  a  lineal  element  of  the  circumference  and  equals  cd8.  e^  is  a 
unit  vector  in  the  radial  direction  and  V  is  the  velocity  vector 


V  *-u  e  +  ua  ea 
r  r  9  0 


(6.28) 


eg  being  the  tangential  unit  vector.  Thus, 

•2it 


Lim 

e-K) 


r 

I  cu  d0 


(6.29) 


where  uf  is  evaluated  at  r«e.  Substituting  for  ur  from  Eq.  (6.25)  into  Eq. 
(6.29),  carrying  out  the  indicated  limit,  and  simplifying,  it  can  be  shown  that 


A 


4ir 


(6.30) 


Introduction  of  above  result  into  Eq.  (6.23)  yields  the  final  expression  for  4> 
for  a  modified  sink  singularity, 
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It  Is  Interesting  to  investigate  Eq.  (6.31)  when  the  bottom  is  flat  and  is  of  unit 


depth,  i.e.,  6-0  in  Eqs  .  (6.11)  and  (6.31).  In  the  limit  6-K),  it  may  be  verified 


Thus,  Eq.  (6.1)  which  represents  velocity  potential  for  constant  bathymetry 


In  susmiary,  the  jet-induced  flow  field  for  a  steady  potential  motion  over  an 


arbitrary  bathymetry  is  governed  by  Eq.  (6.2).  The  velocity  potential  $  in 


this  equation  can  be  shown  to  possess  elementary  singular  solutions  similar  to 


the  classical  solutions  of  Laplace's  equation.  For  the  special  case  of  an 


exponentially  sloping  bottom,  the  elementary  solution  is  given  by  Eq.  (6.31) 


where  the  singularity  is  a  modified  Bessel  function  rather  than  a  logarithm.  Next 


one  should  consider  the  application  of  Eq.  (6.31)  to  determine  the  jet  induced 


flow  field.  This  is  the  subject  of  the  following  section 


6.4  Application  of  Elementary  Solutions  to  Determine  Jet-induced  Circulation 


Patterns  for  Variable  Bathymet 


In  this  section,  an  outline  of  the  possible  use  of  Eq.  (6.31)  to  determine  the 


jet-induced  flow  field  is  given.  First  of  all,  a  physical  interpretation  of  Eq 


(6.31)  is  necessary.  It  is  obvious  from  Eq.  (6.32)  that  the  equipotential  lines 


(i.e.,  lines  of  constant  potential  $)  for  a  sink  are  perfect  circles  and  the 


streamlines  (i.e.,  lines  of  constant  stream  function  ^)  are  rays  emanating  from 


the  origin.  The  streamlines  and  equipotential  lines  are  mutually  orthogonal.  A 


similar  interpretation  of  Eq.  (6.31)  will  be  very  useful  in  understanding  physically 


the  flow  field  characterized  by  the  elementary  solution. 

Future  efforts  will  attempt  to  solve  the  circulation  due  to  a  simple  jet 
(Fig.  3.1)  with  an  exponential  bathymetry.  Following  the  procedure  in  Sec.  3.1, 
the  jet  will  be  represented  by  a  distribution  of  elementary  singularities  in 
Eq.  (6.31),  where  q  must  be  related  to  the  jet  solution  for  an  exponential  bottom. 
Since  the  work  in  Ref.  11  does  not  treat  this  special  case,  jet  solutions  will 
have  to  be  obtained  using  the  general  formulation  presented  therein.  It  is  hoped 
that  these  solutions  will  be  analytical,  but  they  can  be  obtained  numerically, 
if  necessary.  Further,  a  correspondence  can  be  established  between  the  volume 
flux  due  to  the  singularity  and  entrainment  due  to  the  jet  which  will  enable 
expressing  q  in  terms  of  the  jet  parameters.  The  flow  field  can  then  be  deter¬ 
mined  by  integrating  Eq.  (6.31)  with  appropriate  variation  of  q  over  the  entire 
distribution.  The  result  will  be  an  equation  similar  to  Eq.  (3.25).  Velocity 
components  can  be  obtained  by  differentiation  in  Eqs.  (6.24). 

6.5  Effects  of  Unsteadiness  and  Earth's  Rotation  -  Geophysical  Fluid  Dynamics 

Approach 

The  formulation  and  solution  of  jet  induced  circulation  considered  thus  far 
has  ignored  the  effects  of  unsteadiness  of  the  tidal  jet  and  the  Influence  of  the 
Earth's  rotation.  On  a  time  scale  which  is  comparable  to  the  tidal  period,  tem¬ 
poral  effects  definitely  become  important  in  a  tida.'  c.  Moreover,  tidal  jets 
often  extend  to  several  kilometers  offshore  making  the  effects  of  the  Earth's 
rotation  relevant.  In  this  section,  a  preliminary  formulation  is  presented  and 
only  major  steps  in  the  derivation  are  Included.  A  few  possible  avenues  of 
solution  are  also  mentioned  briefly. 
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Figure  6.1  Derivation  of  Geophysical  Fluid  Dynamic  Equations 


Figure  6.1  shows  an  unsteady  water  suface  n(x,y,t)  with  arbitrary  bottom 
variations  represented  by  hg(x,y).  H(x,y,t)  is  the  instantaneous  height  of 
water  column.  The  fluid  motion  is  assumed  to  be  quasi-geos trophic  or  geostrophic 
to  the  lowest  order.  Therefore,  the  first  order  horizontal  velocity  components 
u  and  v  are  independent  of  z.  This  is  equivalent  to  the  shallow  water  approxi¬ 
mation  encountered  earlier  in  Sec.  2.1.  Under  the  assumption  of  quasi-geostrophy, 
the  continuity  equation  can  be  combined  with  the  so-called  kinematic  boundary 
condition  at  the  free  surface  and  the  non-penetration  condition  at  the  bottom 
to  obtain 


This  equation  has  been  derived  in  Ref. (45)  for  the  case  of  a  flat  and  steady 
free  surface  but  the  derivation  can  be  easily  extended  to  account  fora  time- 
varying  free  surface.  Eq.  (6.33)  can  be  interpreted  physically  as  follows. 

For  every  fluid  particle  located  at  an  elevation  'z'  in  the  coordinate  system  of 
Fig.  6.1,  the  ratio  of  its  height  above  the  bottom  to  the  total  height  of  liquid 
column  at  that  location  is  a  constant  as  the  particle  traverses  the  flow  field. 
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Eq.  (6.33)  has  been  derived  from  kinematic  considerations  alone.  Dynamic 
considerations  are  Introduced  through  the  classical  Ertel's  theorem  of  geophysical 
fluid  dynamics.  The  theorem  states  that. for  an  inviscld  constant  density  (and 
hence  trivially  barotropic)  fluid. 


.  „J.  o  (6.3, 

where  C  ■  relative  vorticity  of  fluid,  ft  ■  angular  velocity  of  Earth,  and  A  is 
a  property  that  is  conserved  with  fluid  motion.  From  Eq.  (6.33),  (z-hg)/H  is 
such  a  property  and  it  is  substituted  into  Eq.  (6.34)  in  place  of  A.  Upon 
simplifications,  the  result  is 


D_ 

Dt 


0 


(6.35) 


Now,  due  to  quasi-geos trophy,  the  first  order  motion  is  two  dimensional  and  only 
that  component  of  c  which  is  parallel  to  ft  in  Fig.  6.1  is  significant,  i.e.. 


I  —i  _  3v  9u 

I'l  •  1  *  -17 


(6.36) 


Furthermore,  in  Fig.  6.1, 


2ft  -  2ft  sin  8  -  f 


(6.37) 


where  6  •  latitude. 


Combining  Eqs.  (6.35),  (6.36),  and  (6.37),  one  obtains  the  vorticity 
equation 


(6.38) 
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The  quantity  in  brackets  in  Eq.  (6.38)  is  called  potential  vorticlty 


n 


StL 

H 


(6.39) 


and  thus 


Dll 

Dt 


0 


(6.40) 


Eq.  (6.40)  expresses  the  conservation  of  potential  vorticlty.  It  shows  that 
vorticlty  associated  with  a  fluid  particle  must  increase  when  it  encounters  an 
increase  in  the  instantaneous  depth  H.  Moreover,  even  if  the  particle  possesses 
zero  vorticlty  at  some  instant  of  time,  it  will  acquire  vorticlty  due  to  depth 
changes  solely  because  of  the  presence  of  planetary  vorticlty  f.  Note  that 
these  observations  hold  for  an  inviscid  fluid  and  the  vorticlty  is  not  derived 
from  fluid  viscosity. 

Eq.  (6.38)  can  be  converted  to  a  more  useful  form  in  terms  of  the  stream 
function.  Since  the  lowest  order  motion  is  geostrophic,  which  represents  a 
balance  of  pressure  forces  and  Coriolis  forces,  the  stream  function  is  related 
to  the  pressure  as  follows 


u 


_Ll£  .  _  3(p/ pf ) 
pf  3y  ay 


v  .  J_  .  ajpfr.P. 

pf  3x  3x 


Therefore,  a  stream  function  may  be  defined  as 


(6.41a) 

(6.41b) 


* 


=  JL 
pf 


(6.42) 
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This  is  che  mosc  general  governing  equation  which  allows  free  surface  variations. 
Eq.  (6.44)  is  quite  complex  in  its  present  forn.  Introducing  the  following 
simplifying  assumptions 

(1)  variations  in  H  are  small,  and  therefore,  H  *  Hq  in  Eq.  (6.44) 

(2)  the  relative  vorticity  C  (or  equivalently  V2ty  in  Eq.  (6.43))  is  much  smaller 
compared  to  the  planetary  vorticity,  f,  i.e.,  V2ip«{  in  Eq.  (6.44) 


it  may  be  verified  that 


V2*  +  Ktilhl  .  X  1M  .  V  , 

Vt  3(x,y)  H  3(x,y)  H 


(6.45) 


In  Eq.  (6.45),  the  Jacobian  in  the  third  term  and  the  unsteady  depth  term 
can  be  expressed  in  terms  of  t|»  using  Eq.  (6.42)  and  the  relations 

n(x,y,t)  «  H(x,y,t)  +  h^x.y)  (6.46) 

and 

p(z)  -  patffl  +  pg(n-z)  (6.47) 

The  resulting  equation  in  is 


L 
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(6.48) 


- 


V2h,  +  *>.  +  X 

VC  3(x,y)  Ho  3(x,y) 


fZ 


gH, 


which  can  be  written  more  compactly  as 


+ 


3  0KV2<H-fh_/H  ) 
_ B  O 

3<x,y) 


0 


(6.49) 


The  second  term  in  the  coefficient  of  iji  in  Eq.  (6.49)  owes  its  existence  to  the 
unsteady  free  surface  variations.  In  fact,  this  term  does  not  appear  in  a 
similar  derivation  contained  in  Ref. (45). 


Finally,  through  some  algebraic  manipulations  Eq.  (6.48)  can  be  put  in  the 

form 


D_ 

Dt 


5  f2 

V2’*'  -  — -  4,  + 

L  ®0 


0 


(6.50) 


This  form  is  very  illuminating  in  that  it  shows 

2  c2 


7  i|»  - 


8H 


<l>  + 


‘B 


H 


constant  along  a  path  line 


(6.51) 


Perhaps,  it  is  a  logical  first  step  to  solve  the  steady  flow  problem  including 
the  effects  of  bathymetry  and  Earth's  rotation.  In  this  case, 
fh 

a  2 

V  +  — —  *  constant  along  a  streamline  (6.52) 

o 

Note  that  the  term  in  f2  has  disappeared  since  it  is  a  contribution  due  to 
unsteadiness  as  seen  from  Eq.  (6.49).  The  constant  on  the  RHS  of  Eq.  (6.52) 
varies  from  one  streamline  to  another,  and  therefore,  it  is  a  function  of  ifi 
only.  Thus, 
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(6.53) 


fh 

V2*  +  -jj~  -  F  (ij>) 
o 

Ref.  (46)  provides  a  starting  point  for  the  choice  of  the  function  F.  If  F  is 
known  for  a  physically  meaningful  steady  flow  situation,  Eq.  (6.53)  can  be  solved 
as  a  non-homogeneous  reduced  wave  equation.  In  the  particular  case  of  F  =  0 
and  the  absence  of  bathymetric  variations,  Eq.  (6.53)  reduces  to  Laplace's 
Eq.  (2.12)  which  governs  steady  potential  flows.  Eqs.  (6.53)  and  (6.51)  indi¬ 
cate  that  in  the  absence  of  planetary  rotation  (f  ■  0),  the  bathymetric  effects 
are  also  removed  from  the  problem.  This  is  because  of  assumption  (1)  under  Eq. 
(6.44)  that  the  depth  variations  are  gentle.  In  other  words,  the  effects  of 
small  bottom  variations  are  exemplified  only  when  significant  rotation  is 
present. 

Solution  of  the  general  unsteady  Eq.  (6.49)  is  very  difficult.  References 
(46),  (47),  and  (28)  should  be  very  useful  in  this  regard. 
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7.0  SUMMARY  AND  CONCLUSIONS 


7.1  Summary 

An  analytical  Investigation  of  the  circulation  induced  by  a  tidal  jet 
in  an  ebbing  ocean  was  conducted.  The  jet  was  assumed  to  be  steady  and 
issuing  into  coastal  waters  having  constant  depth  but  arbitrary  bottom 
friction.  The  surrounding  fluid  motion  was  assumed  to  be  steady  and  poten¬ 
tial  (i.e.,  inviscid  and  irrotational) .  The  effects  of  unsteadiness.  Earth's 
rotation,  and  buoyancy  have  been  neglected  in  this  preliminary  study;  however, 
the  Important  effects  of  lateral  mixing  have  been  included.  During  the  final 
stages  of  the  study  reported  herein,  mathematical  formulations  were  developed 
to  describe  steady  potential  flow  with  variable  bathymetry  and  unsteady, 
inviscid,  rotational  flow  due  to  Coriolis  forces. 

To  obtain  the  shelf  circulation  patterns  over  constant  bottom  topography, 
the  jet  was  replaced  by  an  equivalent  distribution  of  sinks  which  are  them¬ 
selves  solutions  of  Laplace's  equation.  The  latter,  of  course,  governs  the 
potential  motion  of  the  jet-induced  circulation.  The  strength  of  each  sink 
in  the  distribution  is  determined  by  equating  the  volume  flow  due  to  entrain¬ 
ment  to  the  volume  flow  demanded  by  the  sink*  Thus,  a  complete  physical 
and  mathematical  connection  was  established  between  the  jet  flow  and  the 
outer  flow.  This  formulation  of  jet-induced  motion  has  been  adopted  from 
Taylor  (38) .  Analytical  expressions  were  derived  for  the  stream  function 
and  velocity  components.  These  quantities  were  evaluated  analytically 
whenever  possible  and  numerically  using  Gauss -Legendre  quadrature.  The 
streamline  patterns  obtained  for  a  simple  jet  (i.e.,  perpendicular  to  the 
coastline)  bear  a  close  resemblence  to  the  circulation  patterns  observed  by 
Llpplsch  (48)  in  a  smoke  tunnel.  The  variation  of  alongshore  current  with 
distance  from  the  jet  was  also  predicted  and  shown  to  be  a  strong  function 
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of  bottom  friction  of  the  jet.  These  numerical  computations  were  performed 
for  various  values  of  bottom  friction.  The  following  four  jet  configurations 
of  practical  interest  were  studied: 

(i)  A  simple  jet  normal  to  the  coastline. 

(ii)  A  jet  with  jetties  of  varying  lengths. 

(ill)  Two  parallel  jets. 

(iv)  A  jet  non-perpendicular  to  the  coastline. 

7.2  Conclusions 

Computational  results  for  the  above  cases  were  generated  to  obtain  the 
jet-induced  circulation  patterns,  the  alongshore  current,  and  the 
onshore-offshore  current.  The  general  conclusions  drawn  from  the  numerical 
results  are: 

Continental  Shelf  Circulation  Patterns: 

(i)  The  jet-induced  circulation  patterns  are  very  sensitive  to  the  bottom 
friction,  f,  and  the  half-width/ depth  ratio,  bQ/ho>  0f  the  inlet.  The  com¬ 
bined  effect  of  these  two  variables  can  be  expressed  in  terms  of  a  single 
parameter  p  -  fbo/8hQ.  For  two  inlets  of  identical  dimensions  and  throat 
velocities,  the  jet  with  less  bottom  friction  (small  f)  Induces  circulation 
over  a  much  larger  area  of  the  continental  shelf  than  a  jet  with  large 
bottom  friction  (large  f).  This  is  because  the  offshore  extent  of  the  jet 
is  reduced  significantly  by  the  Increased  bottom  friction.  Similarly,  for 
two  jets  with  identical  values  of  the  friction  factor  f  and  width  bo>  a  Jet 
with  a  smaller  aspect  ratio  bQ/ho  induces  circulation  over  a  greater  portion  of 
the  continental  shelf  than  a  jet  with  a  large  aspect  ratio. 
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(ii)  The  introduction  of  a  jetty  affects  the  jet-induced  circulation 
patterns  very  characteristically.  The  streamlines  close  to  the  inlet 
display  strong  curvature  toward  the  mouth  of  the  jetty  and  thus  demon¬ 
strate  clearly  the  Influence  of  the  jet-induced  entrainment. 

(iii)  For  the  case  of  two  parallel  jets,  the  induced  flow  field  is  symmetric 
about  the  centerline  between  the  jets  which  divides  the  flow  field  into 

two  regions.  The  waters  in  between  the  jets  are  entrained  by  both  the 
jets.  The  circulation  patterns  in  the  outer  regions,  however,  are  similar 
to  those  for  a  simple  jet. 

(iv)  The  outer  circulation  patterns  associated  with  a  jet  which  is  not 
perpendicular  to  the  coastline  exhibit  asymmetries  (when  compared  to  the 
symmetric  streamline  pattern  for  a  simple  jet)  as  one  would  expect. 

Alongshore  Current: 

(i)  The  alongshore  current  toward  the  inlet  has  a  maximum  value  (about 
3.5  percent  of  the  throat  velocity)  near  the  inlet  and  drops  off  as  one 
moves  away  from  the  inlet  along  the  coastline.  This  drop-off  is  a  strong 
function  of  the  bottom  friction  of  the  jet.  An  increase  in  bottom  friction 
results  in  a  significantly  faster  decay  of  the  alongshore  current. 

(ii)  The  Introduction  of  a  jetty  changes  the  character  of  the  variation  of 
the  alongshore  current  from  that  described  above.  The  current  is  zero  at 
the  inlet,  increases  to  a  maximum  value  at  some  location  along  the  coast, 
and  decreases  to  zero  far  away  from  the  inlet.  The  maximum  value  attained 
depends  upon  the  jetty  length  and  reduces  sharply  with  increased  jetty 
length . 
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(ill)  For  two  identical  parallel  tidal  jets,  it  was  found  that  the  presence 
of  a  second  jet  reduces  the  alongshore  current  (below  the  value  for  a  simple 
jet)  between  the  jets,  but  increases  it  in  the  outer  flow.  This  conclusion 
is  based  on  the  assumption  that  each  jet  entrains  the  same  amount  of  fluid 
as  a  simple  jet  in  spite  of  the  proximity  of  the  other  jet.  Consequently, 
the  entrainment  on  the  inner  boundaries  of  the  jets  drops  and  increases  on 
the  outer  boundaries. 

(iv)  Under  the  same  assumption  regarding  jet  entrainment  as  in  (iii)  above, 
for  a  jet  non-perpendicular  to  the  coastline,  the  alongshore  current  on  the 
acute  angle  side  is  reduced,  whereas  it  is  increased  on  the  obtuse  angle  side. 

Onshore-Offshore  Current: 

The  behavior  of  onshore-offshore  currents  was  investigated  for  a  simple 
jet  and  the  conclusions  are  expected  to  be  similar  for  other  cases,  except 
for  minor  modifications.  For  a  bottom-frictional  jet,  the  current  is  always 
directed  onshore.  As  one  proceeds  from  the  coast  in  the  offshore  direction 
parallel  to  the  jet,  the  magnitude  of  the  onshore  current  increases  from 
zero  to  a  maximum  value  at  some  point  and  drops  off  to  zero  further  offshore. 
The  peak  value  of  the  onshore  current  decreases  as  one  moves  farther  from 
the  jet. 

Steady  Potential  Flow  with  Variable  Bathymetry; 

The  problem  of  steady  potential  flow  with  variable  bathymetry  was  also 
formulated.  It  was  shown  that  arbitrary  bottom  topography  can  be  handled  by 
introducing  a  velocity  potential  rather  than  a  stream  function.  The  sinks 
for  Laplace's  equation  were  replaced  by  special  singularities  whose  form  was 
obtained  for  an  exponentially  sloping  bottom.  It  was  found  that  a  generalized 
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series  solution  can  be  employed  to  determine  the  nature  of  the  singularities 
for  arbitrary  bathymetry.  It  is  expected  that  the  jet  can  be  replaced  by  an 
equivalent  singularity  distribution  and  the  outer  flow  solution  obtained  in  a 
manner  similar  to  that  used  for  the  constant  bathymetry  case. 

Effects  of  Unsteadiness  and  Earth's  Rotation: 

The  effects  of  unsteady  jet-induced  circulation  in  the  presence  of  the 
Earth's  rotation  and  variable  bathymetry  were  formulated  within  the  framework 
of  geophysical  fluid  dynamics.  It  appears  that  the  special  case  of  steady 
and  invlscid  but  rotational  flow  with  arbitrary  bottom  variations  represents 
a  logical  extension  of  the  problem  detailed  in  the  previous  paragraph. 

Recommendations  for  Future  Research  Efforts; 

The  short-term  goals  of  future  research  in  this  area  should  include 
(i)  a  refinement  of  the  algorithm  which  generates  stream  function  and 
streamlines  for  the  constant  depth  case;  (ii)  a  parametric  study  of  the 
effects  of  jetty  length,  distance  between  two  adjacent  jets,  jet  angle  with 
respect  to  the  coastline,  etc.  on  the  circulation  patterns  and  the  alongshore 
current;  (iii)  a  complete  analytical  solution  of  jet-induced  circulation 
for  an  exponential  bathymetry  and  an  extension  to  other  bathymetric 
features  for  which  analytical  solutions  are  possible;  (iv)  a  study  of 
possible  analytical  solutions  to  the  steady,  lnviscld,  and  rotational  jet- 
induced  circulation  with  arbitrary  bathymetry  and  Earth's  rotation. 

The  long-term  goals  of  research  programs  in  this  area  should  be  to 
include  the  effects  of  unsteadiness  in  the  jet  as  well  as  the  induced  flow, 
buoyancy  of  the  jet,  wave  climate,  and  viscosity  in  the  outer  flow.  In  an 
entirely  different  but  parallel  effort,  the  mass  transport  associated  with 
circulation  induced  by  tidal  jets  should  be  investigated. 
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